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We  have  here  considered  a  gas  of  electrons  with  a  positive  background 
interacting  among  themselves  through  a  bare  Coulomb  potential  at  T  -  0°K. 

On  examining  th<*  single-particle  states  of  the  system,  one  discovers  that  the 
bare  Coulomb  potential  is  screened  by  a  dynamic  dielectric  constant.  The 
same  screening  is  shown  to  appear  also  when  the  low-lying  excitations  of  the 
spin-wave  type  involving  particle -hole  pairs  are  examined.  Now,  in  the 

random  phase  approximation  and  in  the  static  limit  where  the  plasma  effects 

/ 

are  neglected  and  for  long  waves,  this  screening  is  the  usual  Thomas-Fermi 
screening.  One  may  thus  start  with  a  Yukawa  potential  which  contains  an 
arbitrary  screening.  This  has  not  only  the  advantage  of  taking  into  account 
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tome  correlations  but  also  includes  the  extreme  long-range  and  the  extreme 

, 

short-range  Interactions  as  special  cases.  The  coupling  strength  is  then 
defined  as  proportional  to  the  ratio  of  the  exchange  to  the  kinetic  energy. 
These  considerations  motivate  the  use  of  the  Yukawa  potential  for  the  present 


problem. 

We  compute  the  ground  state  energy  of  such  a  system  in  the  Hartree- 

Fock  approximation,  but  excluding  for  the  moment  the  possibility  of  spin 

•  * 

density  waves.  This  energy  is  a  function  of  coupling  strength,  .screening 

. 

5  and  relative  magnetization  J  .  Examining  the  absolute  minimum  of  this 
energy  as  a  function  of  magnetization  for  various  coupling  strengths  and 
screening  constants  we  find  some  interesting  features.  For  f  £  0.9  we 
find  the  ground  state  is  either  paramagnetic  or  ferromagnetic  depending  on 
the  value  of  the  coupling  strength;  for  £  >  0.  9  the  unsaturated  ferromagnetic 
states  also  compete.  For  £  <  0.  9  the  intermediate  states  are  relative  maxima 
and,  hence,  do  not  appear.  The  case  where  £  =  0  is  the  Coulomb  problem 
and  was  studied  earlier  by  Bloch  in  1929.  For  £  >  0.  9  and  much  larger 
tending  to  infinity,  one  has  the  extreme  short-range  model  of  Stoner  (1938). 

The  results  obtained  by  us  go  over  into  this  case  smoothly.  Thus,  a  kind  of 
phase  diagram  is  obtained  which  describes  the  various  ground  states  as  the 
screening  and  the  coupling  strengths  are  varied.  It  is  interesting  to  point 
out  that  if  the  Thomas-Fermi  value  is  used  for  the  screening,  it  is  found 
that  the  gas  stays  paramagnetic  for  all  densities. 


The  spin  wave  .excitations  in  the.fystem  are  also  considered  within 

/  . 

the. context  of  the  random  phase  approximation  (RPA).  The  equation  describing 
these  is  solved  for  long  wavelength*  of  the  spin  wave.  The  existence  of  spin 
waves  is  clearly  brought  out  by.  computing  the  coefficient  of  which  is  the 
first  non-zero  term  in  this  expansion.  This  .coefficient  contains  two  terms  with 
opposite  signs.  By  demanding  that  the  spin-wave  frequency  be  positive',  one 
may  define  the  stability  of  the  spin  waves.  The  condition  for  this  coincides  with 
that  obtained  for  the  ferromagnetic  case  by  Herring  both  for  the  Coulomb  and 
for  the  Stoner  gases  when  appropriate  limits  are  taken.  A  consequence  of  this 
criterion  is  that  if  the  ground  state  is  a  stable  ferromagnetic  state  (i.  e. ,  stable 
with  respect  to  individual  particle  excitations)  then  it  is  also  stable  with 

t  * 

respect  to  collective  spin-wave  excitations  in  the  long  wave  limit.  There 
exists  a  maximum  wave  vector  for  the  spin  wave  beyond  which  it  is  unstably 


being  scattered  into  individual  particle  states. 

A  matrix  Green's  function  formalism  is  here  developed  to  treat  this 


problem.  This  formalism  is  a  modification  of  Nambu's  for  superconductivity 
theory.  The  method  is  quite  general  and  is  here  employed  also  to  extend 
formally  the  results  to  electrons  in  periodic  potentials  (Bloch  electrons). 


The  same  method  is  also  used  to  derive  the  equations  for  the  symmetry 


breaking  solutions  of  the  Overhauser  spin  density  wave  type. 

In  an  epilogue,  some  objections  to  the  use  of  Yukawa  potentials  are 
considered.  Also,  the  various  aspects  of  extensions  of  the  problem,  some 
of  which  are  only  formally  developed  in  the  report,  are  outlined. 
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I.  MOTIVATION  AND  STATEMENT  OF  THE  PROBLEM 


J 


1.  General  Introduction 


There  are  essentially  two  basic  models  of  ferromagnetism  in  solids, 
and  all  others  may  be  considered  as  variations  on  these  two  themes.  The 


Heisenberg  [1]  model  is  based  on  the  assumption  that:  (1)  the  electrons  are 
localized  on  atoms;  (2)  only  a  single  electron  configuration  corresponding  to 
one  electron  per  atom  need  be  considered;  and  (3)  the  interatomic  exchange 


effects  can  be  treated  by  introducing  an  interaction  -  2j  S.  *  S,  between  the 

ij—i  — j 

electrons  localized  on  the  sites  i  and  j  .  Here  J\j  is  an  exchange  integral, 
and  Sj  are  the  spin  operators  corresponding  to  the  e’ectrons  at  the  sites 
i  and  j  .  This  model  is  particularly  suited  to  the  case  of  insulators  and  will 


not  be  discussed  here.  For  a  recent  review  of  exchange  in  insulators,  one 


may  refer  to  Anderson  [2].  The  itinerant  model,  on  the  other  hand,  developed 
by  Bloch[3],  Stoner  [4],  and  Slater  [5],  is  based  on  the  competition  between 
the  kinetic  energy  of  the  electrons  in  a  band  and  electron  exchange  in  the 


Hartree-Fock  approximation  (her eafter  referred  to  as  HF).  This  exchange  may 


be  thought  of  as  arising  primarily  from  the  intra-atomic  rather  than  inter¬ 
atomic  exchange.  There  are  two  extreme  models  of  exchange  based  on  long- 


range  Coulomb  interactions  between  the  electx-ons  (which  were  considered  by 
Bloch  [3])  and  on  zero-range  interaction.  The  latter  is  shown  here  to  be 
equivalent  to  a  Weiss  field  and  such  a  model  was  examined  by  Stoner  [4] 
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and  Slater  [!>]  .  Various  modifications  of  the  itinerant  and  the  localised 

o 

schemes  to  fit 'fSfcr  real  solids  may  be  found  in  the  review  articles  by  BroOi>-» 

I 

[6a,  b]  and  Vonsovskii  and  Isyumov  [7]  . 

It  may  be  surmised  that. in  an  itinerant  theory,  the. system  may  become 

ferromagnetic  only  for  a  certain  ratio  of  the  exchange  energy  to  the  kinetic 

energy.  This  ratio  may  be  thought  of  as  a  suitable  strength  parameter,  which 

determines  the  magnetic  behavior.  In  the  next  two  sections  a  brief  summary 

will  be  given  of  the  results  found  in  the  literature  on  the  ground-state  properties 

and  also  the  collective  excitations  of  such  systems. 

» 

2.  Ground-State  Properties 

In  order  to  describe  the  grour.d-state  properties,  the  total  energy  of 
the  system  must  be  computed  and  its  absolute  minimum  as  a  function  of 
relative  magnetization,  ?  ,  must  be  examined.  This  cannot  be  done  exactly. 
The  results  to  be  described  are  within  the  HF  approximation.  This  consists 
in  assuming  that  the  system  behaves  as  if  it  were  composed  of  "quasi” 
particles  which  are  effectively  free  despite  interaction.  Also, this  approximatio 
takes  into  account  only  the  parallel  spin  correlations  between  the  electrons. 

Neglect  of  antiparallel  spin  correlations  tends  to  exaggerate  the  tendency 

v 

towards  fer romagnetismj  that  is,  it  overestimates  the  difference  in  inter¬ 
action  energy  between  parallel  and  antiparallel  spins  ct  the  electrons.  This 
difference  tends  to  be  minimized  when  the  interactions  are  short  range,  and 
therefore,  antiparallel  spin  correlations  can,  to  some  extent,  be  taken  into 
account  phenomenologically  by  introducing  a  short-range  interaction  between 


qua  s’ i  particles.  Most  of  the  calculations  on  the  ferromagnetism  of  an  electron 
-gas  disregard  such  correlations.  It  is  only  recently  that  the  interactions  have  bem 
taken  into  account  more  fully  for  the  ordinary  (nonmagnetic  )  electron  gas.  A 
description  of  the  ground-state  properties  ir.  HF  and  some  recent  attempts  to 
take  account  of  correlations  will  be  given  in  the  next  two  sub-sections.  It  will 
be  assumed  that  the  system  of  electrons  has  a  positive  background  so  as  to 
keep  the  entire. system  electrically  neutral.  All  the  calculations  pertain  to 
T  =  0°K  ,  although  t..  ormalism  will  be  such  that  this  restriction  is  not 


A.  Coulomb  Interaction 


Bloch  [3]  was  the  first  to  compute  the  total  energy  of  an  electron  gas 
in  a  positive  background  at  T  =  0°K  in  the  HF  approximation,  as  a  function  of 
the  relative  magnetization,  J  .  The  coupling  strength  for  this  problem  is  r# 
which  is  the  effective  radius  of  the  electron  in  units  of  the  Bohr  radius  and  is 


inversely  proportional  to  the  cube  root  of  the  density  of  the  system.  He 

rw 

showed  that  for  low  densities, corresponding  to  r  >  5.  4 5 ,  the  gas  becomes 


ferromagnetic.  This  follows  when  the  energy  of  the  ferromagnetic  state 
(  f  =  1  )  is  compared  with  that  of  the  paramagnetic  state  (  £  -  0  )  .  A  similar 
calculation  was  made  by  several  other  authors,  notably  Brillouin  [8], 
Wohlfarth  [9]  ,  Lidiard[10],  Shimuzu[ll],  and,  ve  ry  i  ecently,  Fukuda  [  1 2]  . 
All  obtained  the  same  result  as  Bloch.  Shimuzu  tried  to  take  into  account 


the  electron  correlations  by  including  the  plasma  effects,  but  his  analysis 
was  inconclusive.  Cooper  [13]  extended  the  Gell-mann  and  Brueckner  [14] 


vnc  orainary  (nonmagnetic)  electron  gas  to  the  present  problem. 
He  concluded  that  taking  into  account  more  interaction,  give,  ferromagnetism 
°r  r.  >  metkod  con. lifted  in  taking  into  account  a  certain  infinite 

aet  of  ring  diagram..  The  conclusion  i.  that  the  electron  ga«  with  Coulomb 
interaction,  in  HF  become,  ferromagnetic  at  low  den.itie.  and  paramagnetic 
at  high  den.itie*.  None  of  the  author.  Quoted  abnv»i»m.  _ a 


B.  Short-Range  Interaction 


Stoner  [4]  and  Slater  [5]  modified  the  Somme rf eld  model  of  the  metal 
by  postulating  a  phenomenological  internal  magnetic  field  in  order  to  describe 
the  ferromagnetic  interaction.  The  coupling  strength  in  this  case  is  the  ratio 
of  the  internal  magnetic  field  energy,  K 0',  to  the  Fermi  energy,  eF  ,  K0'/c, 
They  assumed  that  the  internal  Weiss  field  acting  on  each  electron  is  the  sam 
so  that  the  up  and  down  spin  bands  are  displaced  rigidly.  Slater  [5],  however 
took  into  account  the  actual  density  of  states  distribution  instead  of  the  free 
electron  one.  The  second  assumption  in  this  theory  is  that  the  electrons  in 
their  separate  bands  obey  Fermi-Dirac  (FD)  statistics.  Such  a  system  wai 
shown  to  have  a  paramagnetic  (P)  ground  state  for  K0'/e„  <  2/3  an 
unsaturated  ferromagnetic  state  (UF,  0  <  ?  <  1)  for  2/3  <  K0'/cf  <  2"1/3, 
and  a  ferromagnetic  state  (F)  for  K 0'/eF  >  2'1/3  .  Wohlfarth  [15]  has 
recently  reviewed  the  present  status  of  this  model. 

The  outstanding  assumptions,  the  existence  of  a  Weiss  field  and  the 
use  of  FD  statistics,  remained  to  be  justified  in  this  theory.  The  first  was 


.  -rJ-:-, 
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attempted  by.  Wohlfarth  [9]  and  Lidiard  [10]  and  the  second  by  Bell  [16]  . 

Wohlfarth  [9]  and  Lidiard  [10]  showed  that  the  assumption  of  Stoner  that  the 

exchange  energy  is  proportional  to  the  square  of  the  relative  magnetization  is 

, 

not  strictly  true.  More  generally,  it  is  given  by  a  power  series  in  J  ;  the 
coefficients  of  this  series  involve  certain  overlap  integrals.  Earlier, Brooks 
[6b]  had  examined  in  detail  the  ground-state  properties  for  such  a  generalized 
Weiss  field.  Bell  [16]  explicitly,  examined  the  .statistics  of  the. Weiss  field 
and  showed  by  a  direct  calculation  of  the  partition  function  that  the  use  of  FD 
distributions  for  the  up  and  down  bands  can  be  made  consistent.  This  is 


achieved  by  summing,  for  each  magnetic  quantum  number  m  ,  over  the  sets 
of  occupation  numbers  which  are  consistent  with  that  value  of  m  .  Recently, 
Suris  [17]  reformulated  the  Stoner  problem  using  a  Green's  function  method. 

He  found  that  Stoner's  model  is  equivalent  to  an  electron  gas  interacting 
through  a  zero-range  potential  in  the  HF  approximation.  The  FD  distributions 
for  the  up  and  down  spins  follow  as  a  natural  consequence  as  does  also  the 
Weiss  form  of  the  exchange  energy.  Suris  then  reproduced  the  results  derived 
by  Stoner  although  his  method  of  solution  of  the  resulting  integral  equations 


was  more  analytical  than  was  the  case  with  Stoner's  work.  Thus,  this  model 


is  shown  to  exhibit  an  F  state  for  high  densities,  a  P  state  for  low  densities 


and  for  the  intermediate  densities,  UF  states,  in  contrast  to  the  beliavior  of 


the  Coulomb  gas 


All  the  above  remarks  concern  the  symmetry-preserving  solutions  of 
the  HF  equations.  Recently,  Overhauser  [18]  showed  that  there  exist  symmetry¬ 
breaking  solutions  of  the  HF  equations  which  have  lower  energy  than  the  usual 
solutions.  In  particular,  he  showed  that  in  a  Coulomb  gas,  the  P  state  is 
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never  the  ground  state.  It  Is  always  unstable  towards  the  formation  of  static 


spin-density,  wave  (SDW)  states.  This  has  been  confirmed  by  Fukuda  [12]  . 


In  an  electron  gas  with  zero-range  interaction,  according  to  Overhauser,  the 


SDW  states  do  not  occur  as  ground  states  at  all,  for  any  coupling  strength  or 
density.  But  Fukuda  finds  that  these  states  are  more  stable  than  the  P  state 


for  coupling  strengths  for  which  the  UF  is  known  to  be  stable. 


PI! 


:  s 


3.  Individual  Particle  and  Collective  Excitations 

The  low-lying  excitations  in  a  medium  close  to  the  ground  Btate  are 
equally  important  in  the  study  of  any  system.  In  the  Heisenberg  model  of  a 
ferromagnet,  besides  the  spin-aligned  state  constituting  the  ground  state, 
there  exist  spin-wave  states  close  to  the  ground  state.  The  existence  of  spin 
waves  in  both  ferromagnetic  metals  and  insulators  has  been  established 
experimentally  [6a]  .  Originally, the  existence  of  spin  waves  was  proved 
theoretically  by  Bloch[19]  only  for  the  Heisenberg  model.  For  quite  some 
time,  the  absence  of  spin  waves  on  an  itinerant  model  was  thought  to  be  a 
serious  drawback  of  the  itinerant  theory.  Slater  [20],  however,  demonstrated 
that  taking  into  account  interactions  more  fully  would  produce  spin  waves 
on  such  a  model.  This  view  was  later  amplified  by  Herring  and  Kittel  [21] 
and  by  Herring  [22,  23]  .  The  first  paper  gave  largely  phenomenological 
arguments  and  a  more  rigorous  development  appeared  in  the  subsequent 
papers.  In  an  insulator  one  may  describe  the  spin  waves  as  follows.  In  the 
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ground  state,  the  spins  are  all  localized  on  the  lattice  sites  of  the  insulator  and 

o 

are  all  lined  up.  The  first  excited  state  is  such  that  one  of  the  spins  on  one  of 
the  lattice  sites  is  reversed  and  since  the  system  has  translational  symmetry 
this  propagates  as  a  spin  reversal  mode.  Herring  and  Kittel  [21]  gave  arguments 
to  show  that,  despite  the  fact  that  the  electrons  are  wandering  throughout  the 
metal  lattice,  their  . spin  motions  are  correlated  because  of  the  exclusion 
principle.  They  showed  that  an  electron  moving  in  a  ferromagnetic  metal  is 
continually  gyrated  because  the  electrons  in  its  immediate  vicinity  tend  to 
align  the  spin  of  this  electron  with  their  average  spin  moment.  So  they  postulate 
a  phenomenological  expression  for  the  energy  due  to  this  torque  involving  a 
"Bloch  wall  coefficient"  which  is  shown  to  be  related  to  the  spin-wave  energy 
in  the  lowest  order  of  perturbation  theory.  Herring  [22,  23]  treated  this  model 
in  detail  both  for  the  Coulomb  gas  and  for  the  short-range  model.  The  approach 
was  to  treat  the  ferromagnetic  material  as  a  continuous  medium  in  which  the 
three  components  of  spin  density  are  regarded  as  the  amplitudes  of  a  vector 
field  quantized  in  the  way  demanded  by  the  known  commutation  rules  for  the 
spin  components.  They  showed  that  there  exist  low-lying  spin-wave  states  and 
that  these  are  orthogO'.T-l  to  all  the  low-lying  individual  particle  states  of  the 
usual  itinerant  electron  model.  This  same  picture  is  employed  by  all  the 
subsequent  authors  including  the  present  ones. 


ARPA-11 


A.  Coulomb  Gas 


- 


Herring  [22]  coniputed  the  spin-weve  frequency  in  the  long  wavelength 
limit  for  the  ferromagnetic  state,  using  the  model  proposed  earlier  [21]  .  For 
the  Coulomb  gas  he  found  that  the. coefficient  of  q^  (where  q  is  the  magnitude 
of  the  wave  vector  of  the  spin  wave)  consists  of  two  terms.  By  demanding  that 
the  spin-wave  frequency  be  positive,  he  derived  a  stability  condition  which 
showed  that  stable  spin  waves  exist  for  r#  >  5.  485.  This  is  slightly  larger 
than  the  value  derived  by  Bloch  [3]  using  the  ground-state  criterion  (r^  >  5.45). 
Recently,  Fukuda  [12]  rederived  this  result  using  an  equation  of  motion 
technique  in  the  random  phase  approximation  (RPA)  including  exchange,  also 
in  the  limit  of  small  q  .  He  obtained  a  value  of  rg  >  5.  145.  However,  he 
made  a  numerical  error  in  computing  this  value.  The  correct  stability 
criterion  is  rfl  >  5.  344.  This  result  seems  more  plausible  than  Herring's, 

since  it  implies  that  when  r  is  large  enough  to  make  the  ferromagnetic 

s 

state  stable  with  respect  to  individual  particle  excitations  (essentially  the 
Bloch  criterion),  then  the  criterion  for  spin-wave  stability  is  automatically 
satisfied.  In  other  words,  the  Bloch  criterion  is  stronger  than  the  spin-wave 
criterion,  and  as  we  shall  see,  this  result  also  applies  in  the  case  of  the  general 
Yukawa  interaction. 
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B.  Short-Range  Gas 


Herring  [23]  also  computed  the  spin-wave  energy  in  the  long  wave  limit 


for  the  ferromagnetic  state  using  the  short-range  model.  He  computed  the 


coefficient  of  q  as  in  the  Coulomb  case  and  found  that  the  spin-wave  stability 
condition  is  Kfl'/fp  >  2^^/5  =  0.  635.  The  required  coupling  is  smaller  than 
that  required  for  stability  of  the  complete  F  state  (Kfl'/tj.  >  0.794)  or  for 
that  matter,  the  criterion  for  the  UF  state  (Kd'/cp  >  0.  667).  The  same  result 
is  also  derived  by  Fukuda  [12].  This  implies  that  the  spin  waves  are  stable 


whenever  either  UF  or  F  state  is  a  stable  ground  state.  Besides  finding  the 


spin  waves,  Herring  and  Kittel  [21]  had  given  arguments  to  show  that  there 


are  other  excited  states  which  are  of  the  Stoner  type  particle-hole  pair 


scattering  states.  These  were  shown  to  be  orthogonal  to  the  spin-wave  states. 


Hereafter,  we  shall  refer  to  the  spin-wave  states  as  "bound  states"  since  they 


are  states  of  lower  total  spin  polarization  than  the  ground  state  but  lying  lower 


in  energy  than  the  lowest  excited  individual  particle  states  corresponding  to 


the  transfer  of  one  electron  of  one  spin  distribution  to  the  other.  Thompson 
[24]  has  recently  proposed  an  extension  of  the  usual  dete rminantal  method  for 
treating  the  ferromagnetic  metal.  His  work  is  very  close  in  spirit  to  that 
of  Herring.  He  considered  only  the  short-range  interaction  model.  His 
calculation  is  more  general  than  Herring's  in  that  he  computed  the  spin-wave 
dispersion  law  for  the  F  state  for  arbitrary  q  .  Besides  rederiving 
Herring's  expression  for  the  coefficient  of  q^  ,  he  showed  that  the  spin 
waves  become  degenerate  with  the  single  particle  excitations  of  the  Stoner 


t 
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type  above  a  certain  q^^  .  Thio  is  shown  schematically  in  Fig.  la.  This 

degeneracy  actually  results  in  the  breaking  up  of  the  bound  state  into  ths 

scattering  states.  The  spin  wave  decays  by  carrying  an  up  spin  electron  into 

a  spin  down  state,  and  energy  and  momentum  must  be  conserved  in  this 

process.  The  difference  in  energy  of  the  up  and  down  electrons  is  finite  as 

shown  in  Fig.  laf  and  the  momenta  must  he  such  that  the  up  electron  must  be 

inside  the  up  Fermi  sphere  while  the  down  electron  must  be  outside  the  down 

Fermi  sphere.  With  these  restrictions,  for  small  q  ,  this  process  is 

energetically  impossible  while  for  a  certain  q  and  beyond  it  is  indeed 

max 

possible.  This  process  can  take  place  for  all  suitable  momenta  of  the 

electrons  satisfying  these  conditions  and,  hence,  one  gets  a  continuum  of 

scattering  states  as  shown  in  the  figures.  The  bounding  curves  are  for  the 

electrons  near  the  Fermi  sphere.  If  the  calculation  is  followed  carefully,  it 

is  seen  that  a  lifetime  appears  for  the  bound  state  as  soon  as  it  reaches  the 

continuum  of  scattering  states.  The  interesting  feature  is  that  for  ^  <  1  ,  the 

gap  at  q  =  0  ,  2K 9  '  J? ,  decreases  and  so  for  the  UF  states  the  q  is 

max 

bound  to  be  smaller  than  that  for  the  F  state  and  finally  for  f  =  0  ,  the  P 
state,  there  are  no  bound  spin-wave  states.  This  picture  is  essentially  the 
same  for  the  Coulomb  case  also. 

Edwards  [25]  tried  to  extend  and  generalize  Herring'3  work  for  Bloch 
states  but  did  not  report  any  explicit  calculations.  He  derived  an  improved 
version  of  the  general  expressions  for  the  spin-wave  dispersion  in  the  long 
wavelength  limit  using  diagrammatic  techniques  for  the  evaluation  of  matrix 
elements.  He  also  indicated  how  the  coefficients  of  higher  powers  in  q 
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could  bo  calculated.  Kubo,  et  al.  [26]  ,  using  a  Green's  function  technique, 

treated  the  short-range  itinerant  electron  model  in  RPA  including  exchange. 

These  authors  were  interested  in  computing  the  spin^wave  contribution  to  the 

neutron  scattering  for  an  itinerant  electron  ferromagnetic  system.  Such  a 

model  was  shown  to  be  capable  of  accounting  for  the  observed  diffuse  scattering 

In  this  calculation  it  is  necessary  to  compute  the  dynamical  susceptibilities. 

The  spin-wave  dispersion  so  derived  agrees  with  Thompson's  [24]  .  Baym  [27], 

using  a  Boltzmann  equation  technique,  derived  the  same  results  for  the  same 

model.  He  also  gave  arguments  like  Thompson's  for  the  spin-wave  states 

» 

strongly  mixing  with  the  individual  particle  states  above  a  certain  wave  vector. 
It  is  also  of  interest  that,  using  the  Landau  theory  of  Fermi  liquids  as  a  model 
for  the  electrons  in  a  metal,  Abrikosov,  et  al.  [28]  ,  showed  the  possibility  of 
spin  waves  in  a  ferromagnetic  metal.  Antonoff  [29], using  an  equation  of 
motion  method, also  attempted  this  problem;  he  employed  what  he  called  a 
"degenerate  kernel  approximation,"  besides  RPA,and  obtained  spin  waves.  He 
also  computed  the  maximum  wave  vector  qmax  as  was  done  by  Thompson 
and  Baym.  Fukuda  [12]  has  examined  the  problem  of  spin  waves  based  on  the 
Stoner  model  but  he  only  computed  the  spin-wave  frequency  in  the  long  wave¬ 
length  limit  and  in  the  ferromagnetic  state  like  all  others.  Thus  several 
authors  have  computed  the  spin-wave  frequency  in  the  long  wavelength  limit 
and  for  the  F  state  for  the  Stoner  model.  It  must  be  added  that  Kubo,  et  al.  , 
and  Baym  (loc.  cit.  )  found  zero  sound  type  collective  oscillations  to  exist, 
also.  These  results  are  summarized  in  Table  1. 
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TABLE  1  !  Individual  Particle'aTKj- -Collective  Excitations 

(Stoner  Gas) 


Type  of  Disturbance 

density-density 
long,  spin'susc. 
trans.  spin  3usc. 


Collective 

Mode 


Scattering 

States 


U~q^(zero  Cu=E  (k+q) 

a~q) Boun<i)  {(,/(,  !) 

">  / _ N 


,  zero 
lound) 

w~q2  ( spin 
wave) 


~Er(k), 


(zero 

sound) 


J response 


U=  E^(k+q)-E^(k)|no  resonant 

Eff(k)  =  k2/2m-  vNff  (N^  =  k^/6  #2) 

=  radius  of  <r-Fermi  sphere,  v  :  interaction  potential. 


The  density-density  response  function,  longitudinal, and  transverse  spin 
susceptibilities  may  be  thought  of  as  describing,  respectively,  the  spin  singlet 
and  the  spin  triplet  with  projections  0  and  1  of  the  particle-hole  pairs. 
The  scattering  states  for  the  density  response  function  and  the  longitudinal 
spin  susceptibility  are  of  the  form  shown  in  the  schematic  sketch,  Fig.  lb. 
Here  the  zero-sound  mode  is  seen  to  be  quite  close  to  the  scattering  states.' 
The  plasma  state  is  shown  in  this  same  figure,anticipating  the  future 
discussion  of  the  Coulomb  gas.  The  absence  of  an  energy  gap  in  the  excitation 
spectra  must  be  noted. 

It  is  Interesting  to  compare  the  corresponding  results  for  the  P  state 
when  the  interaction  is  taken  to  be  of  zero  range.  Here  only  zero  sound 
appears  in  the  density  response  while  no  resonant  responses  appear  in  the 
susceptibilities  as  might  be  expected.  In  thi3  limit  the  results  are  well  known. 
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(Zero  sound:  Gottfried  and  Pieman  [30]  and  spin  susceptibility:  Wolff  [31]  ). 
The  corresponding  calculations  for  the  Coulomb  gas  are  not  found  in  the 
literature. 

4.  Present  Work 

From  the  brief  summary  of  the  previous  work  in  the  last  two  sections, 
it  is  clear  that  the  itinerant  electron  model  of  ferromagnetism  has  recently 
seen  a  revival  of  interest.  The  conclusions  of  Bloch  [3]  and  subsequent 
authors  on  the  Coulomb  gas  were  slightly  obscured  because  the  absolute 
minimum  of  the  ground  state  was  not  carefully  examined.  For  the  P  state, 
much  work  on  the  many-body  aspects  has  been  accomplished  using  modern 
field-theory  techniques.  For  a  recent  review  of  this  one  may  refer  to  the 
books  of  Pines  (32,  33]  and  Anderson  [34],  Shimuzu  [11]  used  the  Bohm-Pines 
method  to  examine  the  many-body  effects  on  the  criterion  for  ferromagnetism. 
His  results  were  inconclusive  because  he  did  not  examine  the  total  energy  as 
a  function  of  magnetization.  Except  for  the  work  of  Cooper  [  1  3]  ,  the 
extension  of  field  theoretic  methods  to  the  ferromagnetic  problem  has  not  been 
carried  out  and  even  this  work  is  not  without  objections.  One  of  the  main 
objections  to  Cooper's  work,  as  he  himself  realizes,  is  that  the  ring  diagrams 
that  were  summed  are  valid  only  for  high  density  (r  <  1)  and  the  second  major 
objection  is  that  he  only  compared  as  Bloch  did,  the  F  and  P  energies. 

It  has  often  been  suggested  in  the  literature,  notably  by  Lidiard  [10] 
and  Wohlfarth  [35]  ,  that  some  account  of  correlations  can  be  included  by 
using  a  Yukawa  potential  for  the  interaction.  This  contains  a  screening 
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parameter  and,  hence,  include*  both  the  extreme  long-range  Coulomb  cate 
and  the  extreme  short-range  .Stoner  case  as  limits.  The  suggestion  actually 
came  from  Landsberg  [36]  who  tried  to  explain  the  soft  x  '.ay  emission  bands 
of  sodium  on  a  Sommerfeld  model  of  the  metal.  Very  recently,  Robinson  et.  al 

[37] ,  considered  the  justification,  within  the  context  of  RPA,  of  the  use  of 
such  a  Yukav/a  potential  for  nonmagnetic  problems.  The. screening  can  bs 
shown  to  come  about  when  particle  interactions  are  taken  into  account.  Pines 

[38]  had  earlier  shown  that  taking  into  account  the  electron-electron  inter¬ 
action  brings  about  a  screened  interaction  between  electrons,  of  range  k  ” ^ 

(in  his  notation).  He  showed,  even  though  this  has  a  very  different  structure, 
that  it  resembles  the  Yukawa  potential  quite  closely  up  to  distances  comparable 
to  or  less  than  kc  .  All  these  arguments  pertain  to  the  nonmagnetic  case. 
Whatever  be  the  consideration  for  the  use  of  Yukawa  potential  in  an  inter¬ 
acting  system,  it  seems  fruitful  to  re-examine  the  problem  afresh  as  the 
behavior  of  the  Coulomb  and  the  Stoner  gases  are  completely  different.  The 
Yukawa  potential  may  give  an  insight  into  the  effect  of  range  on  the  criterion 

of  ferromagnetism. 

The  use  of  the  Yukawa  potential  may  be  justified  for  the  magnetic 
problem  in  much  the  same  way  as  was  done  by  Robinson,  et  al.  [37]  .  If^ 
in  an  interacting  electron  gas,  the  behavior  of  a  single  electron  is  examined, 
it  is  seen  that  even  though  one  starts  with  a  bare  Coulomb  interaction  the 
exchange  is  screened  in  general  by  the  dynamic  dielectric  constant.  In  the 
RPA,  and  if  plasma  effects  are  neglected  by  assuming  the  dielectric  constant 
to  be  static,  and  working  in  the  long  wavelength  limit,  it  is  found  that  the 
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interaction  potential  ha.  now  the  appearance  of  a  Yukawa  potential  with  the 
screening  parameter  given  by  the  familiar  Thoma.-Fermi  screening.  This 
same  screening  can  be  shown  to  appear  when  the  collective  excitations  of  the 
spin-wave  type  are  examined.  In  this  way,  the  use  of  the  Yukawa  potential 
may  be  justified.  This  will  be  amplified  in  later  chapters. 

The  Green's  function  method  is  employed  here  as  it  is  found  to  be 
powerful  enough  to  include  all  the  results  of  previous  authors  and  also  to 
extend  them,  as  will  be  indicated  in  subsequent  sections  .  Thus  a  unification 
of  all  the  results  both  as  to  the  nature  of  the  interactions  and  as  to  the  use  of 
a  unified  and  elegant  technique  for  attacking  the  present  problem  is  here 
attempted.  All  the  previous  authors  concerned  themselves  with  plane  waves 
while  the  present  work  indicates  extensions  to  Bloch  waves  also. 

Besides  making  a  more  complete  study  of  the  ground  state  of  a 
polarized  gas,  the  few  non-rigorous  attempts  of  Wohlfarth  [15],  Lidiard  [10], 
and  Bell  [  1  6]  to  justify  the  use  of  a  Weiss  field  and  FD  distributions  will  be 
here  given  a  more  satisfactory  justification.  This  consists  in  showing  that 
a  quasi  -particle  picture  is  valid  in  metals  as  was  done  by  Luttingcr  [39] 
for  the  ordinary  interacting  electron  gas,  whereas  the  authors  quoted  above 
could  only  show  the  consistency  of  such  assumptions.  They  could  prove  that 
the  exchange  energy  is  in  general  a  power  series  in  Here  it  is  seen 

to  be  a  complicated  function  of  ?2  ,  and  so  the  convergence  becomes  slower 
as  ?  approaches  unity.  Bell  showed  how  the  partition  function  could  be 
correctly  evaluated  in  the  presence  of  a  Weiss  field  and,  hence,  also  established 
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the  validity  of  the  FD  distributions  for  the  particle.  These  fall  out 


automatically  in  the  present  approach. 

On  the  collective  aspects  of  the  sytem,  most  of  the  authors  confined 
themselves  to  the  short-range  Stoner  model  and  the  F  state.  For  completenes 
the  Coulomb  gas  must  also  be  studied  in  the  same  spirit.  This  problem  is 
much  harder.  In  the  present  work  the  spin  waves  are  studied  in  detail  in  the 
long  wave  limit  and  a  formal  power  series  expansion  Is  also  given  for  the 
frequency  of  the  spin  wave.  Herring's  [22,23]  method,  though  ingenious,  is 
valid  only  in  the  long  wave  limit.  He  used  a  perturbation  method  for  the  HF 
equations  and  from  the  agreement  with  the  present  work,  this  may  be  viewed 
as  RPA  in  a  different  form.  Thompson's  [24]  dcterminantal  method  is  very 
cumbersome.  The  method  of  Edwards  [25]  gives  a  technique  for  evaluating 
the  coefficients  of  various  powers  of  q  and,  hence,  is  also  quite  cumbersome. 
The  method  of  Kubo,  ct  al.  [26],  uses  a  factorization  of  the  Green's  functions 
which  is  equivalent  to  RPA,  and  this  is  here  rederived.  Baym's  [27] 

Boltzmann  equation  technique  is  equivalent  to  RPA  right  from  the  start  and 
it  is  hard  to  see  the  complete  equation  before  the  approximations  are  made. 
Antonoff's  [29]  method  is  similar  to  that  of  Kubo,  et  al.  j  his  use  of  the 
"degenerate  kernel"  approximation  may  be  shown  to  be  equivalent  to  the 
localization  of  the  spins,  with  a  Heisenberg-type  interaction.  This  invalidates 
the  purpose  of  a  fully  itinerant  theory.  The  same  objection  is  also  applicable 
to  the  work  of  Shimuzu  [11],  who  explicitly  as s umes  in  an  ad  hoc  fashion 
a  Heisenberg-type  interaction  between  the  electrons. 
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Thus,  it  is  seen  that  there  is  only  a  sketchy  account  of  the  polarised 
gas  using  field  theory  techniques  in  the  litersnEttreb^There  is  thus  scope  for 
further  generalization  of  the  problem  in  that  one  could  include  the  results 
found  in  the  literature  as  special  cases  of  the  Yukawa  potential. 


A.  Ground-State  Properties 

The  ground-state  energy  of  a  system  of  electrons  in  a  positive  back¬ 
ground  in  HF  approximation  at  T  =  0°K  is  here  computed  as  a  function  of 
?  (0  <?  <  1  )  when  the  electrons  are  interacting  through  a  Yukawa  potential. 

The  problem  is  studied  in  detail  as  a  function  of  screening  (£),  magnetization 
(  ?  )  and  density  (  arg  )  .  To  determine  the  ground  state,  the  absolute 
minimum  of  the  total  energy  as  a  function  of  ?  for  various  arg  and  1-  is 
studied.  It  is  found  that  for  f  much  smaller  than  the  diameter  of  the  ferro- 

l 

magnetic  Fermi  sphere,  the  system  behaves  essentially  like  the  long-range 

I 

Coulomb  system.  It  is  found  in  this  case  that  the  F  state  is  the  ground  state 

i 

for  low  densities  and  the  P  state  for  high  densities.  The  UF  states  can 
never  lie  lowest  as  they  are  relative  maxima.  This  indicates  why  the 
comparison  of  just  the  energies  of  F  and  P  states  in  the  Coulomb  case  is 

I 

sufficient.  For  ?  much  larger  than  the  diameter  of  the  ferromagnetic 
Fermi  sphere,  the  ga3  behaves  like  the  short-range  Stoner  gas  with  all  the 
three  states  of  magnetization  having  the  possibility  of  being  the  ground  state. 

This  difference  in  the  behavior  comes  about  because  of  the  different  forms 

! 

the  exchange  energy  takes  in  these  two  limits.  In  this  connection,  it  may 
be  mentioned  that  Brooks  [6b]  ,  using  a  generalized  Weiss  field  model 


consisting  of  a  power  series  in  f  of  the  exchange  energy  (  cf.  Wohlfarfch  [15]) 

had.  shown  that  the  UF  states  cannot  be  ground  states  if  the  ratio  of  the 

4  2 

coefficient  of  J  to  that  of  J  exceeds  a. certain  critical  value.  This 
observation  had  escaped  the  notice  of  Wohlfarth  (Loc.  cit. )  who  did  not  examine 
the  consequences  of  the  generalized  Weiss  field.  The  parameter  corresponding 
to  Brooks’s  here  is. the  screening.  It  is  clear  from  the  present  analysis  that 
as  £  increases  the  criterion  for  the  F  state  changes  abruptly  beyond  a 


certain  critical  value. 


These  results  are  summarized  in  Fig.  2a  t  b. 


Here  a  plot  of  the  coupling  strength  (In  ( arg)  in  Fig.  2a  and  K0'/cp  in  Fig.  2b) 
against  $  is  presented.  This  clearly  shows  that  for  £,<,0.9  the  nature 
of  the  ground  state  changes.  A  more  detailed  description  will  be  given  in 
Section  II.  The  possibility  of  SDW  being  ground  states  has  not  been  examined 
here  for  the  Yukawa  gas.  From  the  work  of  Overhauser  it  is  known  that  the 
ground  state  is  an  SDW  state  for  densities  at  which  the  present  theory  shows 
the  P  state  to  be  lowest. 


*fC 


S  SCREENI 


B.  Collective  States 


effective  masses*  This  is  shown  in  Table  me  transverse 
bility  also  fails  to  show  spin  waves  when  exchange  is  omitted; 


Individual  Particle  and  Collective  Excitations 


(The  Coulomb  Gas) 


Type  of  Disturbance  Collective 


density -density 


no  resonant 
response 


trans.  spin  susc 


radius  of  v  Fermi  sphere;  for  small  k  , 
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equation  describing  it  when  exchange  is  included  can  be  solved  in  the  long 


wave  limit.  It  is  for  this  reason  that  is  Table  2  the  spin  waves  are  sfiSw*^ 
within  brackets.  This  was  also  done  by  Fukuda  [12j  .  The  method  cf  Fukuda 
1  j  used  here  and  it  is  found  that  he  made  an  error  in  the  computation  of  the 
coefficient  of  q2  ;  he  found  in  the  F  state,  ~  (q2/*«  )  (I  -  5. 145/r  ) 

whereas  it  is  shown  here  that  to  the  same  order  of  approximation  as  Fukuda 

*^SW  ~  t q  Z2"1 )  H  -  5.  344/r#  )  .  This  result  must  be  compared  with  that 
derived  by  Herring  [22]  :  «sw  ~  (q2/2m)  (1-5.  485/r9). 

For  the  sake  of  completeness  the  results  for  the  P  state  are  also 


given  in  the  Table.  Here  only  the  plasma  mode  exists  when  exchange  is 
neglected  and  there  are  no  other  excitations.  The  way  the  scattering  states 
and  the  bound  states  are  disposed  is  similar  to  that  given  earlier  in  the 


section  <?n  the  short-range  gas. 


*■  1 

In  the  Yukawa  case,  the  most  interesting  collective  mode  is  the  spin 


flip  mode,  after  including  exchange  contributions.  For  any  J  ,  and  for 


small  q  ,  by  using  the  same  technique  as  Fukuda,  the  coefficient  of  q2  can 


be  determined.  For  the  F  state,  for  £  =  0  (Coulomb)  and  for  $  -  oo 


(short  range),  this  goes  over  to  the  results  obtained  earlier.  In  the  short- 


range  case,  a  general  dispersion  law  for  spin  waves  of  all  q  ,  and  any 
magnetization  f  ,  (at  T  =  0°K)  is  obtained.  The  coefficient  of  q2  for  the 


UF  state  is  here  derived.  For  the  F  state,  this  goes  over  to  the  result  of 


Herring  [23]  .  Also,  the  spin  waves  merge  into  the  scattering  states, which 
are  of  Stoner  type, at  about  3/4  the  F  state  Fermi  momentum  and  the  spin. 
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wav*  at  fcbi*  value  of  the  wave  vector  is  about  l/lfe  of  the  F  state 

Fermi  energy.  .Similar  estimates  were  made  by  Thompson  for  the  F  state 
and  by  Baym  for  the  UF  states  (toe.  cit. ).  In  the  UF  state,  the  spin  waves 

are  .stable  for  *>/»/«  — *  I).  *?l/  *1  — 1 - - 


seen  that 


p  g  g  -i 

are  stable  for  Kfl'/tp  >(l/5)  |  — ^■■■*.  i.\  *{L,  ,j  whereas  the  Stoner 

condition  is  K0'/eF  =  (l/2)  01— &  ~  **  "  ^  ^  j  .  Thus  it  is  seen  that 

the  spin  wave  stability  condition  is  always  weaker  than  the  Stoner  condition. 

Also  the  q  for  the  UF  states  becomes  smaller  as  t  decreases  as 

max 

described  earlier. 

The  formalism  given  here  can  be  modified  to  treat  the  Overha  user 
(18]  problem  of  the  SDW  states.  We  have  rederived  the  Overhatiser  integral 
equations  by  our  method  in  a  straightforward  manner.  The  extension  to  the 
Bloch  states  is  also  indicated  here.  The  collective  excitations  of  the  plasma 
type  for  the  Bloch  electron  system  in  the  unpolarized  case  were  first  derived 
by  Ehrenreich  and  Cohen  [40]  who  used  an  equation  of  motion  method.  This 
is  here  generalized  to  the  polarized  case.  The  other  response  functions, 
namely  the  longitudinal  and  transverse  spin  susceptibilities  are  also  computed 
for  the  Bloch  scheme.  The  possibility  of  the  existence  of  spin  waves  in  the 
Bloch  scheme  could  only  be  indicated  in  certain  crude  approximations  due  to 
the  very  complicated  nature  of  the  equations  describing  them. 

In  summary,  the  Yukawa  gas  behaves  like  a  Coulomb  gas  for  £  <  0.9 
with  no  UF  states  and  for  £  >  0.9,  like  the  Stoner  gas,  as  displayed  in 


Fig..  2a,  b.  Al.o  .pin  wave,  are  pre.ent  on  .uch  a  model,  and  the  coefficient 
of  q  for  the  .pin-wave  frequency  i*  computed.  The  re.ult.  go  ever  into 
those  derived  by  Herring  for  the  Coulomb  and  Stoner  ca.e.  when  appropriate 
limit,  are  taken  (with  some  modification,  for  the  Coulomb  cafe).  It  is  shown 
that  whenever  the  HF  ground  atate  is  an  F  .tate,  it  al.o  exhibit,  .table  .pin- 
wave  excitation..  In  the  Stoner  ga.,  the  .pin  wave,  for  the  F  .tate  are 
stable  when  the  ground  .tate  I.  the  UF  .tate.  Al.o,  it  is  shown  that  there 
exist,  a  q _  for  the  snin  W2V«C  at  tf  ■  mu  iv\ a  m.va  »..!  >L  •  y  ..  >  ■  « 


ARPA-1  1  II.  SINGLE  PARTICLE  EXCITATIONS  AND  TK« 


1.  Introduction 


After  establishing  the  notation  briefly,  a  discussion  of  the  single 

*  *  _  -  •  *  '  • 

particle  .states  of  the  system  under  consideration  will  be  given.  The  notation 

will  be  for  the  most  part  that  of  Martin  and  Schwinger  [41]  or  equiyalently 
Kadanoff  and  Baym  [42]  ;  the  latter  will  be  henceforth  referred  to  as  KB  . 
For  the  sake  of  completeness  a  brief  summary  of  Chapters  1 ,  2  ,  and  5  of 
KB  will  now  bs  given. 

The  quantum  mechanics  of  a  system  of  identical  particles  are  best 
described  in  terms  of  second  quantized  operators,  namely  the  creation 
operator  ij/  ^  (  r,  t,  )  and  the  annihilation  operator  ^  {  r,  t,  )  (when  acting 

to  the  right).  Here  stands  for  the  spin  orientation  and  (  r  j  tj  )  represent; 
the  space-time  point,  which  will  be  represented  in  the  future  by  1  .  The 
dynamics  of  the  system  is  described  by 


+  j  2,  J  J  d  rd  r'dt'  (rt)  (r't')V(r-r';  t-t')^,  (r't'ty/^frt) 

(2.  i.  1) 

in  units  where  jfi  =  1.  Here  the  first  term  represents  the  kinetic  energy, 
the  second  term,  the  single  particle  potential,  and  the  last  term  represents 
the  two-particle  Interaction  potential  between  the  particles.  The  interaction 
potential  is  taken  to  be  instantaneous  so  that 
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.VifLA 

. +  . 


Wa  will  henceforth  write  t-*t  +  0;  G  {  rt ;  r  t*  >  =  G  (  rr‘ ;  0  )  *  G  (11*  )  . 
V{  r-r1  )  la  normally  the  Coulomb  potential  and  it  ie  left  unspecified  for  the 
present  so  as  to  preserve  maximum  generality.  Since  we  are  considering  an 
electron  system,  the  Pauli  exclusion  principle  is  written  in  the  form  of  a  set 
of  anti-commutation  rules  at  equal  times  for  these  operators:  "  " 

[tf/  ^(rt),^  t(r‘t))  =0;  [^r(rt),  (r't)]+  *  «rrl«(3)(r-r'>  (2.1.3) 

9  v  ..  t.  •  .  .  |  „  i  L'.'.t  -  >•-  •  *r 

where  (A,  Bj+  =  AB  +  BA  .  The  equation  of  motion  for  any  operator  X(t) 

in  the  Heisenberg  representation  is  ...  ,:•/  . 

i  ±|Jli  =  [X(t),H)_  .  ,  (2.1.4) 

where  (A,  B]  =  AB  -  BA  .  In  order  to  study  the  thermodynamics  of  the 
system  the  grand  canonical  ensemble  average  is  defined  by 

<  X>  .  Tr  {.-«•  tH-“N)  x}  /Tr  {.'P  '"‘""’l  -  «*•  » 

Here  Tr  stands  for  the  trace  on  the  states  of  the  system;  H  is  the  total 
Hamiltonian  and  N  is  the  number  operator  £  \J/  ^  (  rt  )  ^  (  rt )  d3r  ,  p  is  the 
chemical  potential  and  {5  is  l/KT  with  K,  the  Boltzmann  constant  and 
T  the  temperature  on  the  absolute  scale. 

We  now  define  single  particle  correlation  functions 


(2.1.4) 


G  >  ,  (11')  =  -f  <^(1)  > 

<r«r  i  a  « 

g < ,  (n')=  4-  <^0to') > 

aa'  io  o 


These  are  respectively  defined  for  tj  >  tj,  and  tj  <  tp  and  for  real  times. 
These  functions  define  the  propagation  of  disturbances  in  which  a  single 
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particle  is  added  to  the  extern  at  some  space-time  point  1  and  removed  at 
*■  in  the  first  case  while  the  second  defines  the  opposite  order  of  operations 
The  single  particle  Green's  function  may  be  defined  as 


where  T  is  the  Wick  time  ordering  symbol  for  Fermi 


It  follows  that 


All  the  above  definitions  hold  for  real  time  domains.  Using  the  equation  of 
motion  for  xjsj  1),  the  equation  for  0^(12)  may  be  constructed.  Before 
discussing  this,  a  few  general  remarks  about  G  will  now  be  made  which 
follow  by  virtue  of  its  definition.  Now  formally,  since  Tr  (ABC)  =  Tr  (BCA) 
Tr  (CAB),  one  has  a  relation  between  G >  and  G<  at  the  bonn^rinc 


imaginary  time  domain 


Or,  if  the  domain  of  definition  of  G  is  extended  to  imaginary  times 
0  <  it  <  p  ,  one  may  write  formally 


along,  with 


These  may  be  written  equivalently.  In  the  other  time  variable  also.  Since  the 
system  under  consideration  has  time  translational  invariance,  the  functions 
G>,  G<  ,  G  are  functions  of  the  difference  (tj  -tp)  only.  In  view  of  this, 
(2.  1.9)  may  be  written  as 


Now  the  Fourier  transforms  in  time  may  be  introduced 


i  in  the  first,  and  -i  in  the  second  relation).  Using 
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rSfitr 


so  that  formally  on«  gets 

'  ~  *■: 

C/»,<rlr2i  ">  *  G  *  <r  *  j*>. 

1  c  *2 

introduce  the. spectral  weight  function  A 

rlr2  1  2 

r,»2(rlr2i  ">  *  C\.z  (rlr2i  “>  -  Gr1»2,r-l'**  "> . 


(2.  1.  12) 


We  now 


(2.  1.  13) 


.Clearly  then 

~  > 


(rlr2iU)  =  (1  -  nr(u))A  r  (*i*:2j») 

«  i  z 


} 


(2.  1.  14) 


G<rltr2(rl  rZ>  “)  =  -  nF(u)  A,r1a2(rlr2i  u) 

where  nf  (o)  =  l/(exp  (w-/i)p+l). 

One  may  write  a  Fourier  series  representation  of  G  (1  2) 
of  (2.  1.10)  (Chapter  3  of  KB,  Pp  19,  20) 


in  view 


V2<r.i Tz' ti  tz) = ( ztf)  ^ 


v  (odd 
integers) 


-u  (t,-<2) 

*  \> 


(0  5  itj  ,  i  t2  <  P)  (2.  1.  15) 

where  =  (tt  v/-ip)  +  /i  .  Since  the  condition  (2.  1.  10)  holds  for  both  the 
time  variables,  the  Fourier  coefficient  i3  given  by 

+iz  (t.-tO 

G(rl  <r2^ ^ =  J  C  G«r1a2(rlr2»  tl't2)dtl  • 

This  must,  however,  be  independent  of  t2  and  letting  t?  -  0+  ,  so  that 


and,  hence,  aftor  some  manipulation, 

%  A  (r  r-;u) 

r-  i  »  _  f  d«  *1*2  1  2 

ff,1r2^rlr2*  *y  J  hi  u-*v  ' 

Thus,  the  Fourier  coefficient  is  just  an  analytic  function 


(2.  1.16) 


G  (r.r,  j  z )  =  f 
tr  j<r^  1  2*  '  J  hi 


U)  -  z 


evaluated  at  z  =  z^,  provided  A  has  no  other  singularities.  A  careful 
analysis  of  this  analytic  continuation  is  0iven  in  KB  .  It  is  obvious  that 


A_  _  (r.r  ,u)*Lt  [G  (r  r7;uiit)  -G  (r.  r?;  u-ie)]  ,  (2.1.17) 

*1*2  1  2  e-0  *1*2  1  2  *1*2  1  2 


This  immediately  leads  to  the  sum  rule 


J  ZiT  Acr1cr2(rlr2iu)=^0I  WHC,'G«r1«r2(rlr2;  U'iC)] 

(2.  1.  18) 


=  «  a  fi(3) (rrr2) 

(T  j  cr  *  « 
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becaus* 


1  l%^rlrV  0+>  -  Gc 


,ir2(rlr2i  °‘>1  '  #r  T  «(3|(r,-r2) 


which  is  the  canonical  commutation  rule  at  equal  times.  The  advanced  and  the 
retarded  Green’s  functions  could  equally  well  have  been  introduced 


=  1  ^or  t  >  0  and  0  for  t  <  0  .  These  are  related  to  the  causal 

Green's  function  defined  above  in  a  very  simple  way,  being  equal  to  the 

proper  analytic  continuation  of  G^.  ^  (r^r^;  z)  in  their  respective  domains 

1  2 

of  definition. 

To  simplify  the  notation  of  Green's  functions  G  (12)  may  be 

ala2 

written  in  the  form  of  a  matrix 


G  (  1 2  )  = 


Gj  j  ( 12)  Gj  j  (12) 


Gi  f  (12) 


|  (12) 


(2.  1.  19) 


This  form  resembles  that  used  by  Nambu  [43]  in  his  formulation  of  the 


superconductivity  theory  in  which  the  off-diagonal  Green's  functions  are  the 


anomalous  Green's  functions.  Let  (  tq  ,  Tj  ,  r2  ,  r3  )  be  the  unit  and  the  three 
Pauli  matrices 
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{o  I)'  T‘  =  f  o);  T**P  o )  ^  T’*(l  !,)  <2' 


1.20) 


Then  the  total  number  of  particle#  is  by  definition  n(l)~<^  (1)^|  (l)+^(l)^| 
In  view  of  (2.  1 . 7,  8) 

n{l)  =  -i  [tr  r0G  (11+>J  .  (2.  1.21a) 

.Similarly, the  total  spin  moment  is 

a(l)  =  -i  [tr  r  G(ll+)  )  (2.1.21b) 

with 

<^(1)  =  -i  [tr  Tl  G(ll+)]  =  <ty\  (1)^  (1)+  ^(l)tf'f  (1)>  etc. 

Here  tr  denotes  the  trace  on  the  spin  indices.  The  spin  magnetic  moment  is 
given  by  -jr  £  »  where  is  the  Bohr  magneton. 

The  equation  of  motion  satisfied  by  G(ll')  (henceforth,  no  special 
symbol  will  indicate  matrices  are  being  used,  unless  otherwise  stated)  can 
now  be  constructed,  3inee  the  Heisenberg  equation  for  (1)  is 


‘Vv,,= 


-zsr  *  V<J) 


t  i 1 >  +  } ^ 1  -  3)  J  ( D 


Here,  in  the  third  term  on  the  right-hand  side,  the  times  of  3  and  1  are  the 
same.  Then, 

2 

(12)  - 


i  JL_  +  _! —  _  v  ( l  > 

8tj  2m  1  ' 


alr2 


-X I  SsVa-vl  <T(^j3 >  -  «v 

"T 


(1)>. 


6(4)(l-2) 
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5  <4*(1  -2)  =  $ ( £j  -  r  2)  .  5  {tj  -t2).  The  laat  term  in  the  left-hand  side  of  this 
equation  can  be  recast  in  terms  of  G  and  its  variational  derivative  which  will 
be  defined  presently.  (For  details  refer  KB  Chapter  5).  An  extra  source 
term 

H'  =  J  J  d3r  Ur(rt)  ^(rt)  ^r(rt)  (2.1.22) 

r 

is  added  to  the  Hamiltonian  H  ,  and  this  is  made  to  vanish  at  the  end  of  all 
calculations.  The  operators  are  still  kept  in  the  Heisenberg  representation 
and  the  traces  are  also  taken  over  the  grand  canonical  ensemble  pertaining 
to  the  total  Hamiltonian.  We  shall  not  go  into  the  redefinitions  of  G  in  the 
presence  of  U  and  we  shall  not  even  indicate  by  any  symbol  that  G  is 
evaluated  in  the  presence  of  U  unless  otherwise  stated.  In  view  of  this,  the 
following  well-known  relationship  is  obtained:  (which  may  also  be  taken  as 
a  definition) 

5  1  1 

1  TUT2J"  <x(1)>  =  <t.(^T(2)^(2)X(  1)> -  <^*(2)^(2)  ><X(1)>(2.  1. 23) 

since  U  generates  if/^  .  This  follows  directly  by  developing  <X(1)>  as  a 
power  series  in  U  and  examining  the  term  which  is  linear  in  U  .  Now 

tJ,  <T «'i<3> kjj> k  (i> *<2> )> 

oj  3  12 


can  be  rewritten  by  using  (2.  1.  23)  as  follows: 
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y  i - 2- - J-  <  T  (<b  (1)  d  (2)  >  > 

4  <su.J5;  ‘r'l  '  1  * 

£  <T(^a)^(3)^i(i)42(2>>>- £  <4(3»^»<T(viin.^,2,, 


rj 

or  equivalently 


or 


i 


y  — 2 — _r,(  i2)  =y  4-<Tii//y(T)v!'(DC' 

^  6U  (J;  »,-«.♦>  f  ‘  *T  **  r» 

cr  j  (Tj  '3  1  ft 


^  (2)1  > 


(l)fcT(2» 


+  [itr  t  G(T7+)jG(12) 


4-  y  <T(i//  *  (?)tf/  ($*  (1)^(2))  > 

i  £  r<rj  ffj  cr2 


-y 


5  G  ( 1  2) 


5U  (3;t~=t  ) 

aj  <rj  3  1 


— .  -  [itr  T  G(33+  )]  G  ( 1  2) 
•  %  o 


Thus 


-  v<1»)  To  +  iTo  I  T.co^nj  cu2> 

(4r)  ul(1)  +(-P)  G«‘2)-iI 

- — *  fT'tr 


6G  (12) 


=  T  6 (4>  (1-2)  . 
o 


r 


l 
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We  shall  henceforth  use  .r  *  1  or  assume  its  presence.  An  inverse  matrix 

°  \  ' 


Green's  function  may  be  formally  introduced: 
j*  G_1(12)  G(l3)d42~=  8<4*(l-3) 


JgOIJG"1^)*4!*  8<4Ni-3) 


}■ 


mm  * 


i 


r,, 

! 


(2.  1.24) 


By  taking  the  variational  derivative  of  the  second  of  these,  and  multiplying 

the  resulting  expression  from  the  right  by  a  suitable  G  ,  and  using  the  first 

> 

of  the  relations  in  (2.  1.  24)  we  obtain 


AG 

6 


GUiL.  =  .  C  d44  d45  G (1?)  gG~  l*5),.  G(?2)  .  (2.  1.  25) 

U  (5)  J  6  U  (5) 


We  define  a  new  function,  (12;  3)  which  is  often  called  the  vertex  part 

2) 


o 

-1, 


T  (12;  3)  =  Y  -6G--ii 

°  ^  6U  (3) 


*3  °3 


(2.  1. 26) 


Then 


-  (^r1)ut<‘>  -(-T2)  U|(utij‘d4J^i.5)(lrG(3T+)] 


L 


G(  12) 


+  i  y  d4I  d44  d4?'^l-D  G(14)ro(  4  5;  T)  G(?2)  =  6{4)(l-2)  (2.  1.  27) 


ft 

L. 


j,  . ■  v.  ■!v-b,«s  « .w..y  • 
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Multiplying  this  on  the  right  by  G-1(21')  and  integrating  over  2  and  using 
(2. 1.  24)  we  get 

G  Vl')-  6(4>{1- 

+  i  J  d4I^(l-J)[trG(J5+)J  C(4)(l-l.) 

+  i  Jd^d4?^!-!)  G(l?)  T  (?!•;  J>  (2.1.28) 

The  solution  of  this  equation  must  obey  the  condition  (2.  1.  10)  when  suitably 
extended  into  the  imaginary  time  domain. 

The  new  function  Tq  may  now  be  determined  by  writing  an  equation 

for  it  using  its  definition  and  the  expression  (2.  1.  28)  for  G"1  .  This  equation 

will  now  involve  another  new  function  61^/6  U  and  so  on.  Thus  a  chain  of 

equations  is  obtained.  The  equation  for  T  is 

o 

r0(12;  3)  =  -  6(  *0-2)  6*4*(1  -3)  -  ij* d4Td42d43^/(l-T)tr  [G(l2)ro(23;  3)G(3t+)] 

6(4)(l-2) 

-  ij  d4Td42d43d44t/(l-T)  G(13)ro(34;  3)G(42)T  (?2;1) 


+  ii  d4ld42^(l-l)  G(  1  2)  Tol  (T2;  T3) 


(2.1.  29) 


where 


ARPA-11 


r  .  <  2  2;  T3 )  ■ 


The  various  terms  in  (2.  1.  28)  may  be  interpreted  as  follows.  The  last  two 

t't 

terms  are  due  to  the  two-particle  interactions  whereas  the  remaining  terms 
arise  from  the  a  ingle-particle  parts  in  the  Hamiltonian  and  the  external  source. 
The  first  of  the  last  two  terms  is  the  Hartree  self-energy,  which,  in  the  absence 
of  external  density  disturbances,  cancels  precisely  with  the  positive  back¬ 
ground.  The  last  term  contains  the  contributions  from  all  kinds  of  scattering 
beginning  with  the  exchange  processes.  (The  Hartree  term  may  also  be 
thought  of  as  arising  from  direct  scattering  processes).  In  view  of  this 
interpretation  of  the  terms  in  G  *  ,  the  terms  on  the  right-hand  side  of 
(2.  1.  29)  for  ro  may  be  similarly  interpreted.  The  first  term  is  due  to  the 
direct  external  perturbation,  the  second  ter  n  arises  from  the  density 
fluctuations  stemming  from  the  Hartree  seif  energy,  the  third  term  includes 
all  kinds  of  complicated  scattering  processes  and  the  fourth  contains  various 
forms  of  corrections  arising  from  changes  in  the  vertices  due  to  interactions. 
This  last  term  is  of  second  order  in  the  interaction  potential.  The  equation 
for  r*o  is  seen  to  be  nonlinear  and  inhomogeneous.  Because  of  this,  as  a 
crude  first  approximation,  its  solution  may  be  taken  to  be  just  the 
inhomogeneous  term.  This,  in  (2.  1.28)  gives  the  well-known  Hartree-Fock 
approximation. 

Before  discussing  the  various  forms  of  these  equations,  a  schematic 
outline  of  the  method  of  solution  of  (2.  1.  28)  will  be  given.  Let  U^.  =  0  in 
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(2.^28)  and  the  last  two  term*  in  it  be  denoted  by 


2j  (n,)  =  iJ  (l-T)  [tr  G{TT*)J  $*4,(l-l») 

*i  j'd‘4rd42y(l-7)G(U)ro(21*i  T)  . 

Here  the  Hartree  energy  i.  retained  just  for  the  sake  of  completene...  From 

its  structure,  it  may  be  written  in  terms  of  (  1  .  T,  v  -r  \ 

•1*  2*T3* 

^Tdl')  =  a(  11')  +  jr  JM1I') 

containing  four  unknown  functions.  When  this  is  substituted  back  in  (2.  1. 28), 


G'l(n,)  ’  (l  §17- +  4r  -  v<i>)«,4>(i-i'>+ 


fi.01’). 


Hence  G  is  expressed  in  terms  of  the  unknown  functions  a  ,  p  .  Putting 
this  back  in  £  ,  which  contains  G  ,  fully  determines  the  unknown  function. 

So  far  no  mention  has  been  made  of  the  geometric  structure  of  the 
medium.  Three  cases  wilt  be  considered.  In  the  first  two,  the  single: 
particle  potential  V(r)  will  be  taken  to  be  identically  sero,  so  that  plane 
waves  are  the  solutions  for  the  nonjnteracting  system.  Here  there  are 
two  possibilities.  One  is  the  translationally  invariant  case,  while  the  second 
is  the  symmetry  breaking  solution  of  the  Overha  user  [18]  type.  In  the  third 

case  V(r)  is  taken  to  be  a  periodic  potential  such  that  the  noninteracting 

system  is  now  described  by  Bloch  waves.  Here  again  there  are  two  case. 
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corresponding  to  the  ones  for  the  plane  waves  but  only  the  symmetry  preserving 
solutions  will  be  here  considered.  (We  need  to  know  the  nature  of  the  non- 
interacting  system,  for  the  solutions  of  this  are  used  as  a  complete  set  of 

•  T  7 

states  when  the  system  is  interacting  in  much  the  same  way  the  normal  modes 
of  vibrations  of  a  lattice  are  employed  In  treating  the  enharmonic  vibrations). 


2.  Spatially  Uniform  Solutions 


If  the  system  is  taken  to  be  spatially  uniform,  then  the  Green's  functions 
depend  only  on  the  space-time  coordinate  differences.  In  view  of  this  one  may 
take  Fourier  transform*  in  all  the  variables;  the  Fourier  transform  variable 
corresponding  to  the  space  coordinates  has  the  significance  of  wave  number 
for  single  particle  states.  Let,  therefore. 


We  will  often  make  use  of  the  four-dimensional  notation:  p.  1  =  (j>.  r^ -p^tj)  ; 
P  =  (R>P()^»  po  is  a  frequency  or  energy  variable  (which  will  be  equivalently 
referred  to  as  o  ).  Single  G  is  a  2x2  matrix,  it  can  be  written  in  terms 


of  unit  and  Pauli  matrices 


Correspondingly,  g(|>u)  and  cr(j>w)  are  defined  similarly  to  (2.  2.  1).  Taking 
t.  ,  =  t.  -  0,  we  have  from  (2.  i .  6) 
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With  this  in  {?,.  2.  1)  we  may  define  (c.  f.  equations  (2.  I.  21a,  b)) 

n(p)  =  j*  gl^u)  eiw°  =  total  number  of  particle*  of 

|  momentum  £  (2.2.3) 


(*  du  jmji 

ff(p)  =  j  r(f  u)  *  ww  =  apin  density  of  momentum  j> 


The  simplest  approximation  is  to  take  the  first  term  in  (2.  1.  29)  a*  the 


solution  for  T 


The  expression  (2.  1.  28)  for  G”  after  taking  U  =  0  (also  V(  r)  =  0  in  the 
present  calculation)  takes  the  well-known  Hartree-Fock  form:  (here  the 
Hartree  self-energy  term  is  dropped  as  it  is  cancelled  by  the  positive  back- 


- 


(2.  2.  4) 


ground) 


ghf 


(‘  Tq-  +  zk~) 


Taking  Fourier  transforms, 

ghf  1(j?u)  =  u~p2/2rn"i^'  — ^3— —  2Aj>-]>)  ghf  <E“>  eiU° 

=  u-£2/2m-iC  -  E  ”  JP  >  nr  ^  W  )  AHF  ("pl3 )  e1U°  (2.2.6) 


jr{£«)  =  s(j>)  [g+  (j>«)  -  g_  (j>u)l  J  (■ 

On.  integration  over  the  u  plan*  a*  specified  by.  (2.  2.  3)  the  following 

•  . 

expressions  are  obtained: 

n(jp)  =  (nF(w+{p^))  4-  Hj.  («_(£))  ]  "] 


with  iij.(x)  =  Fermi  function.  The  above  results  were  derived  by 
Suris  [17]  also.  Thus,  defining  n(j>)  and  «r(j>)  as  in  (2.2.3)  and  recalcu- 
lating  them  in  the  HF  approximation  gives  us  the  equations  (2.  2.  11)  to 
determine  them.  Note  that  u+  contain  (n,  a)  .  These  equations  are  thus  • 
nonlinear  integral  equations  for  n ,  r  of  the  Hammerstein  type.  Suris 
(loc.  cit.  )  used  the  theory  of  these  equations  to  derive  the  usual  Stoner  results 
Discussion  of  a  few  aspects  of  these  results  may  be  in  order  here. 

The  up  and  down  spin  states  are  pushed  apart  and  the  energy  difference  is 


This  is  proportional  to  the  total  magnetization  if  one  assumed  'if(p) 
independent  of  £  ,  or  *1/  to  be  short  ranged,  since  I  tr(p)d3p/(2*) 
the  total  magnetization  of  the  system.  From  (2.  2.  11)  then,  one  has 


+  nF(£2/2m  4VN4  (2.  2.  12b) 

These  are  precisely  the  expressions  assumed  by'Stoner[4]  inhls  collective  electron 
theory  of  ferromagnetism.  Thiis  in  the  case  of  zero-range  interaction,  and 
in  the  HF  approximation,  Stoner's  assumptions  are  fully  justified.  The 
Stoner's  constant  K0'  is  seen  to  be  related  to  HT  through  the  relationship 
;  2?  =  (2K9-A0  .  (2. 2.  13) 


These  results  were  derived  by  Suris  [17]  and  also  anticipated  to  some  extent 

by  Thompson  [24]  .  They  were  included  here  for  the  sake  of  completeness 
• 

and  to  provide  a  frame-work  for  generalization  to  the  Bloch  case  and  to  the 
symmetry  breaking  solutions  of  the  SDW  type.  Before  deriving  the  Overhauser 
equations  with  the  present  formalism,  the  relations  (2.  2.  3')  will  be  recast 
in  a  form  that  corresponds  to  the  Overhauser  equations  in  the  special  case 
that  the  symmetry  breaking  parameter  Q  becomes  zero.  To  see  this,  let 
us  define 


J.  (-R.)  ’=  p  -  q )  jr  (  q)  d3q/(2j  )' 


(2.  2.  14) 


(this  |s  actually  related  to  Overha  user's  g(p)  )  ,  Multiplying  the  second  of 
(2.  Z.  11)  by  *l/(p -q)  and  integrating  with  respect  to  j»  ,.we  arrive  at 


Thus 


3.  Symmetry  Breaking  Solutions-  of  the  Overhauser  Type 


Besides  the  usual  plane  wave  solutions  of  the  HF  equations;  Overhauser 
[18]  pointed  out  that  there  are  symmetry  breaking  solutions  which  are  also 
possible  ground  states  of  the  system.  Overhauser  associates  a  plane  wave 
of  the  form  exp  ( ijk  •  r_)  with  the  up  spin  and  exp  (i  ( jk  +  Q  )  •  r)  with  the. 


down  spin.  Thus  the  vector  Q  is  a  measure  of  the  inhomogeneity,  or  in 
other  words,  a  measure  of  the  breaking  of  translational  symmetry.  The 
unrestricted  HF  solutions  of  this  type  are  alternative  candidates  for  the 


ground  state  of  the  system  equally  acceptable  as  the  conventional  ones. 
Overhauser  first  introduces  such  a  term  into  the  HF  equations  and  then  finds 
conditions  for  consistency  very  similar  in  structure  to  (2.  2.  15).  His  g(k) 
is  seen  to  be  a  convolution  of  the  interaction  potential  with  the  intrinsic 


spin  density.  To  derive  Overhauser's  expressions,  the  following 
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redcfiattion.  of  the  Fouri.r  triform,  of  a*  Grce„<,  function.  »„  „.cc...ry 


KWf, 


*  ' *  J  -  ouw.-t.j]  p 

f  d^'  Ifr’+S)- (£,-!,,)  S  PHI 

J  (2»)3  *  G| 

g,,(u')-  f  d3-  ‘k'  f.-Kk'tO)-. 

11  j  uT)1-  GM*i‘r*i.» 

~  •  f  d3K*  i(k!+Q)r.-ik**rlf 

if  ,  =  J  (2*)"1"  *  "  “  G|  f Ck*;4!-*!,)  (2.3.1) 


Then  the  HF  equation  (2.  2.  5)  take,  the  form 


-SJ  f  (is) 


gf|(k) 


u  -  ej  (k+  Q) 


(2.  3.2J 


/here 


-f  <J?  )  =  k2/2m  +  i  f  2/(  k  -  k'  )  G  «  f  (  k'  u'  )  eiu'°+  -  ~ 


^|  (  k  +  Q  )  =  (k  +  Q)2/2m  +  i  jV<  k  -  k'  )Gj  i%'  u’ )  eiu' °+  —Z 


gt  l(-}  =  if  Z'(k-k'  )  Gi^k*  w»)  eiw'0+  —Z 


g|  f  (k)  =  i  J  ?/(k-k')  G.  |VU')  eiu'0+ 


d^k'du* 


_  xX 


Hence,inverting  we  obtain 


(k+<3) 


G(ku) 


where 


g,  (k<j)  =  [  l/(u-u,  (k)  )  ) 


From  (2.  3.  5)  one  has  the  following  easily  verifiable  relationship* 


+  g|  }(JS)  g{  f(ii) 


(k  +  Q)  ) 


(u+(  k)-  ej(  k  +  Q)) 


(u.{k)  -  e.  (k)  ) 


■  ■  y  Y 
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(U+(k)-W_{k))»-2fl - 


€j(  jt)  -  Cj(HQ) 


(2.  3.  ?c| 


(«+{k)  -  c|  (k)  )  («+(k)  -  Cj(k  +  Q)  )  *  f|  j(k)  g}  j  (k)  (2.  3.  7d) 

(w+(k)  -  c|  (k)  )2  +  g|  j(k)  g|  f  * 

.  •  t  ’  •  I 

*  (u+{k)  -  (k)  )  (o+(k)  -  cj(k)  ) 

►  (2.  3.  7«> 

(wjk)  -C|(k)  )2  +  gjj  (Mgjf(k)  = 

~  (u+(k)  -  u_  (k)  )  (u_  (k)  -  e  j{k)  ) 

as  also 

(u+(k)-cj(k  +  Q))2+g}  j(k)gj  j(k)=  (u+(k)  -  w_(k))  (u+(k)  -cj(k  +  Q))  ' 

(u  (k)-  ej(k  +  Q))2+gf  j(k)gj  j(k)  =  (u+(k)-u_(k))  (u,  (k)-Cj  (k  +  Q))  . 

(2.  3.  7f) 

From  the  coefficients  of  g+(ku)  and  g  (ku)  in  (2.3.6)  it  is  reasonable 


to  set 


CO3^0. 


•  za 
sin  v. 


u+(k)  -  C|(k  +  Q) 
u+  (k)  -  u  (  k) 


u  (k)- ci  (k  +  Q) 


(2.3.6) 


whose  sum  is  clearly  unity. 
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'  -~-r  - 
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From  (2.  3.  7e,  f)  the  choice  of  the  branches  u,  may  be  removed  and 

* 


(2.  3.  7ey  f)  may  be  written  as 


* 

-vT 

v 


(w-  cj{k))*  +  g|  |(k)  gj  j(k)  =  (u+<k)  -  o_(k»  («  -ej  <k)) 


(u-ej  (k  +  Q))  *gf  j{k)fijj  (k)  *  (w+(k>  -  v.  {k»  («“  <J  (k  +  Q 


rl!*»  "I 

M2.  3.  7g) 

(w-ei  (k  +  Q))J 


to? 


' 


Using  these  results  for  ( 1c )  -  u_(k))  in  (2.  3.8), 


2  (u+(k)-Cj{k+Q))  (w+(k)-e|(k)) 


CO.  9k  . 


(u+(k) -ej{  k))  +  gf  |(k)  gj  j(k) 


and  in  view  of  (2.  3.  7d) 


cos  ^9. 


gf  |<k)  8J  }<k> 


(  u+ (k)- cj  (k))  +g||(k)g|j(k) 


We  could  define  quite  generally 


2„ 

cos  a. 


g|  |  (k)  gj  |  (k) 

(u-  €|  (  k  ))^  +  gj  |(k)  g|  |(  k  ) 


(2.  3.9) 


with  the  specification  of  the  branch  and  the  definitions  of  Q,  the  same  as 


given  by  Overhauser.  Moreover,  from  (2.  3.8)  and  (2.  3.  7d) 


a.  ^  r  - 


■r 


* 


U\71 
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coo  »k  sin  * 
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’  '  - 

■ 


(u+(k>  r^_(k))2 


v 

) 


_ gf|<fc)  8|  f  (k) _ 

/edlO-eitk  +  QjN2  1 

l - Tf - y  +gf  |i!i)g|  f(]k)J  (2.3. 


We  may  now  write  G(ku)  as 


2CW 


G  (ku)  = 


[co82dkg+(ku)  +  sin20kg(ku)]  coo  0k«in0k  /fLLz^(g+(ku).g  (ku))\ 

-  “  ~N  g|  |(k) 


coafl.sin  0,7— - (g.(ku)-g  (ku»  [sin20,  g  (ku)+cos20,  g  (kw»/ 

-  -  g|  |(k)  +“  '  -  i 


(2.  3.  11) 


These  may  be  substituted  back  in  (2.  3.  3)  to  give 
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I 

I 


■  .  .  .. 


,  i  ,  v.  £  8  ■  0ft 

2  f*  d  k1  _ 

€|(k)  =  k  /2m- \  — ^  ^(k-k')(co.20k,nF(M+(h'))  +  .in\,nF{W  (k‘))]  ^ 


C|  (k  +  Q) 


(k  +  Q)2  r  d3k* 


“  J  ^-^^{k-k'^in^^^Oc’))  +  cos^.nytujk'))] 


C  d3-  / gn  (k1) 

!n  (-,xj  ^^(k-k-jy-^r  ^^•inVnF{w+fri>>‘nF^-^) 

r  d3~  ✓  /in 

(k.y  co^k.«in0k,[nFK(k'))"nF(w_(k 


git(-)aj  w)  co#0k'8in0k*fnFK(?i,))"nFK(ii,))]J 

(2.  3.  12) 

where  the  consistent  sign  of  the  square  root  iri  (2.  3.  10)  must  be  used.  If 
in  these  we  take,  as  was  done  by  Overhauser,  g|j  =  gj|  =  g  ,  and  assume 
only  u^(k')  is  occupied,  then  we  recover  all  of  Overhauser's  results 
including  the  definitions  of  cosS^  .  Thus  a  generalization  of  the  Overhauser 
results  is  here  achieved.  Moreover,  if  we  take  Q  =  0,  we  see  immediately 


U+  (  }  =  2rrT 


l  f  d  — ' 

~  ~T  )  - - T  2/(k  *  k1)  n  (k1  )  +  |s(k)| 

J  (Zl)  ~  ~  - 


with  s  (k)  defined  as  in  (2.  2.  14),  which  are  the  same  as  those  in  (2.  2.8). 
Also  from  (2.  3.  10) 


2-  2  gfl  (k)  gl  ♦  (k) 

cos  0  sin  ti  =  — Li - - 

-  5  4  |  s  (k)  |2 


cose,  sin 9.  y8H  8H 


s 


which  is  the  Tj  component  of  the  equation  (2.  2.  15),etc.  Thus,  the  Overhaussr 
results  are  the  generalizations  of  the  usual  Stoner  conditions  written  in  a 
slightly  different  form.  Incidentally,  (2.  3.  12)  generalizes  the  Overhaussr 
results  for  finite  temperatures  as  well  as  for  finite  magnetization. 


4.  Homogeneous  Solutions  Involving  Bloch  Electrons 


Here  a  brief  description  of  the  extension  of  the  above  results  in  the 

conventional  case  (Q  =  0)  will  be  given  for  Bloch  electrons.  Let  b  ,  (  1  ) 

fk 

be  the  stationary  Bloch  ^unction  satisfying  the  SchrcSdinger  equation 


where  f  is  the  band  indc<  and  k  ,  the  reduced  wave  vector  of  the  electron 
The  Green's  function  in  the  Bloch  k  space  can  now  be  constructed  (  a 
continuum  of  k  is  here  used  for  convenience)  thus 


Here  I  BZ  implies  that  the  sum  on  k  is  over  those  in  the  first  Brillouin 
zone.  This  is  now  substituted  in  the  HF  equation  (2.2.  5),  and  neglecting 


2/00 


where  U<p(  r)  is  the  periodic  part  of  the  Bloch  function,  b 
consistency  conditions  (2.  2.  15}  now  take  the  form 


This  represents  a  system  of  equations  in  which  the  band  polarizations  are 
coupled.  The  neglect  of  Umklapp  processes  is  justified  if  the  bands  are 
broad  and  the  local  field  corrections  neglected  (Adler  [44]  ).  It  must  also 
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be  pointed  out  that  in  deriving  (2.  4.  3)  one  neglects  the  contribution*  from  other 
determinantal  states  (or,  the  configuration  interactions  are  neglected).  These 
appear  in  the  present  formalism  since  the  second  quantisation  technique  includes 
all  the  determinants  that  can  be  constructed  out  of  b^O)*  Tho  properties  of 
the  Bloch  functions  are  such  that  the  matrix  elements  of  the  two-particle  inter¬ 
action  potential  has  a  restriction  only  on  its  k  vector  parts  and  not  on  the 
band  indices.  This  implies  that  there  are  contributions  to  the  matrix  elements 
from  all  the  four  band  indices  and  only  three  wave  vectors. 


5.  Discussion  of  the  Ground-State  Properties 

In  the  last  few  sections  the  single -particle  Green's  functions  for  a 
polarized  interacting  electron  gas  are  studied  in  the  HF  scheme.  This  shows 
that  the  singLe  -particle  states  in  HF  split  into  two,  corresponding  to  the  two 
directions  of  quantization  of  the  electrons,  parallel  and  antiparallel  to  the 
internal  polarization.  In  the  present  section,  the  ground-state  energy  of  the 
system  is  computed  using  the  Green's  functions  derived  so  far..  It  can  be. 
shown  that  the  expectation  value  of  the  Hamiltonian  (2.  1.  1)  with  U  =  0,  can 


be  written  as 

<  H  > 

-  i 

Lt 

I 

H 

ft 

o 

1 

oS 

1  -  1' 

‘rv 

+  0* 

l  Btv) 

4. 

V  ^ 

(  1  - 

yj 

G(ir; 

1  3 

)>  d  r^  volume  of  the  system) 

T 

y  2m 

2m  ) 

(2.  5.  1) 

This  is  related  to  the  free  energy,  F  ,  of  the  system  through  the  relations!-!}. 


(pF)  *  <h >/a 


It  is,  however,  the  ground-state  energy  when  T  =  0°K.  We  shall  discuss  the 
homogeneous  and  the  Ovcrhauscr  cases  separately  (the  latter  is  included  only 
for  the  sake  of  completeness). 

(a)  In  the  spatially  uniform  case,  using  the  Fourier  transforms  as  in 
(2.  2.  1)  for  G(ll'),  and  performing  the  indicated  operations  in  (2.  5.  1)  one 
finally  arrives  at 


Here  we  have  used  a  four-dimensional  notation  as  explained  earlier.  From 
(2.  1. 27)  (with  U  =  0,  V  =  0)  it  is  easy  to  verify  after  defining  the  Fourier 
transform  of  T  to  be 


£<et  us  call  the  last  two  terras  in  brackets  in  the  left 
^which  is  a  2x2  matrix.  Then  formally  one  has 


G(p)  *  - ‘ - 

PQ  -  P  2/2m  +  £  <£  P0) 

Putting  this  in  (2.  5.  3)  and  after  some  manipulation,  one  gets  finally 


This  is  quite  general  and  involves  no  approximations.  The  first  term  is  the 
kinetic  energy,  the  second  is  the  Hartreo  energy  and  the  third  term  contains 
energy  due  to  all  the  interactions  save  the  direct  Hartree  term.  The 
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<H>. 


d3a  a2 

=  J  7r~i isr  I«r^#»rto.iitf)]. 


d3£d3  r  v, 

--rll  ~g-"/^£-a)  J[®F(«+(a))+»r(».(a»3  lnFKte,,+nF{u-(E>>J 
+  i  <e>  *  i<s>  InF  («+<a))  -  nF<«_ta»J  [nF(»+(E»  -  nF(uj[>))]J  (2*  5. 7) 


Suris  [17]  also  derived  an  expression  for  the  total  energy  in  this  form. 

So  far  no  restriction  has  been  made  on  the  form  of  the  interaction 

potential,  2/<q)  .  It  will  be  shown  in  the  next  chapter  that  T  { p  j  p-q)  is 

^  • 

proportional  to  the  inverse  of  a  propagating  dielectric  constant.  In  the  RPA, 
this  dielectric  function  can  be  explicitly  computed.  In  the  static  limit  and  for 
a  long  wavelength,  this  gives  rise  to  a  Thomas  -  Fe  rmi  screening.  The  same 
screening  reappears  when  the  collective  excitations  are  studied.  Thus,  for 
the  purpose  of  the  present  investigation  we  shall  use  a  Yukawa  interaction 
potential  with  arbitrary  screening  parameter,  and  study  the  nature  of  the 
solutions  as  a  function  of  this  screening  (c.  f.  Robinson  et.al.,[37]  ). 


We  take  the  Yukawa  potential  in  the  form 


r  )  =  e  (exp  -  ?  kf  r  )/  r  , 


where  k^  is  the  paramagnetic  Fermi  momentum,  £  is  a  dimensionless 
parameter  and  e  is  the  square  of  the  electronic  charge.  We  also  assume 
that 

nF  (u+  (q))  =  rj+  (kF|  -  q) 
nF  ^q^  =  *4  ^kFJ  '  q^ 
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kj.j  and  k^j  are  given  by 

=  U  +  r)1/3k  ' 


M  =n  +  n^kF-| 
<f|  Mi  -  n1/3kFJ 


where  {  is  the  relative  magnetization  and 

(kpj  ♦  k^j  )/6f  2  =  N  =  total  number  =  k^?  /it  2 

i«  the  positive  unit  step  function.  Furthermore,  the  G's  are  assumed 
to  be  diagonal,  equivalent  to  having  the  internal  polarization  in  the  z 
direction.  The  single-particle  energies  in  the  HF  approximation  then  have 
the  following  form  (c.  f.  (2.  2.  8))  at  T  =  0°K  : 

U+(R)  =JP2/2m  -  e2/*p 


(2.5.10) 

Evaluating  <H>/flQ  under  the  same  approximations,  one  ha3  for  the  total 
energy  per  electron  in  the  HF  approximation 
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with 


,5/3 


Ct  «  f  >  =  4  (i.)  -  ^  ( JL)(1  +  f)4/3^(5/2(i  ♦  C)1/3  ) 


C|  <  ?)  =  -5- 


3  (1  -  g)5/3  •  2 


<^r>  *ET(-4rHl-f)4/3^(f/2(l-»l/3)  J 

(2.  5.12) 


ar#*  ?F  1  ar# 


where  we  have  introduced  ar#  =  me2AF,  a  =  (4/9*  )1^3  =  0.  521,  and 
j/ ( z  )  =  -  2/3  z2  -  8/3  z  tan-*  (l/x)  + 

+  2/3  z2  (z2  4-  3)  fn  (1  +  1/c2)}.  (2.  5.  13) 

^U)  has  the  following  features.  It  is  unity  at  z  =  0  and  monotonically 
decreases  to  zero  as  z  increases,  and  for  z»  1,  ^  (  z  )  ~l/9  z2  ,  which 
is  obtained  by  expanding  3~(  z  )  in  a  power  series  of  l/z.  A  schematic 
diagram  of  its  behavior  is  given  in  Fig.  i. 

The  coupling  strengths  must  be  defined  properly  when  the  two  extreme 
limits  5  =  0  and  5  =  00  are  considered.  These  are  displayed  in  Table  3. 
From  this  table  it  is  clear  that  one  may  infer  the  Coulomb  behavior  for 
C  1  and  the  short-range  behavior  for  5  1.  The  precise  value  of  5 

for  which  this  transition  takes  place  will  be  gi'  in  later  in  this  section.  Also, 
the  coupling  strengths  in  the  two  limits  are  seen  to  be  related  through  the  last 
expression  in  the  table.  Let  us  redefine  the  energy  per  particle  in  terms  of 
the  lermi  energy  so  that  we  only  need  to  examine  a  dimensionless  quantity: 


The  absolute  minimum  of  this  as  a  function  of  the  magnetization  for  various 
values  of  ^  now  be  examined. 


TABLE  3  The  Yukawa  Potential  and  the  Various  Limits 

°0(t)  =  e2{exp  -  ?kp  r)/r  ;  V  (q)  =  4*  e2/(q2  +  £2  kp2  ) 
kj.  :  Fermi  momentum 

For  ?  «  1  ,  M,i  =  4jr  e2/q2  Coulomb  potential 

We  define  the  coupling  strength  as  g^  =  ar( 

(r  in  Bohr  units)  (a  =  (4/97T)1^3  =  0.521) 

For  ?  »  1 ,  °U(  q)  =  47re2/£2  k^.2  =  a  constant 

(Stoner)  zero  range  potential 

We  define  the  coupling  strength  g  =  mTP  k.-./2f  2 

s  r 

=  (2ars)Ae2 

In  terms  of  the  Stoner  coupling  strength  from  (2.  2.  13) 
gg  =  3/2  (K0'/e  F )  ,  ~  Fermi  energy  (5  =  0) 

Thus 

K d'/eF  =  4/3  (ar8/ir?2  )  . 


As  6(C)  *•  defined  on  &  dosed  interval,  the  end  points  must  be  discussed 

» 

separately.  For  the  interior  points  0  <  £  <  1  one  has  the  usual  derivative 

. 

condition.: 


The  first  condition  gives 


xG(5/2x)-yG{?/2y) 

{0  <  f  <  1) 


where 


G(z)  ,  like  j(t.)  ,  is  monotonically  decreasing  with  increasing  z  ,  taking 
the  value  unity  for  z  =  0,  and  going  to  zero  like  l/6z^  for  large  z  .  A  plot 
of  this  is  given  in  Fig.  3.  In  the  above  x  -  (1  +  f)*^  and  y  =  (1-  {)  ^  .  We 
have  to  examine  the  sign  of  8^C/8£^  which  is  given  by  (2.  5.  18)  to  ascertain 
the  existence  of  a  minimum  in  this  interval: 


min 


2 


(2.  5.  18) 


N  .Vj 


Before  discussing  these  la  detail,  a  remark  concerning  (2.-S.  16)  m*/ 
-be  in  order.  It  must  be  noted  that  the  definitions  .(2.  5.  8)  and  {2.  5.  9) 

7  '  j 

the  relations  (2.  2. 11)  a  tautology.  But  if  kj|  and  k^i  are  identified  with 
the  common  chemical  potential  for  the  two  spins  by  the  relationship 

- 

IS  *  «.  (k-*  )  =  u  ( kc_ >  ) 


then  (2.  5.  16)  follows.  Ilera  u+  areas  in  (2.  5. 10).  This  shows  that  t! 
definitions  (2.  5.  9)  are  consistent.  This  will  be  touched  upon  again  late 
The  following  trivial  lemma  is  of  great  use  in  deciding  the  sign 


&2e/d  c2 . 


Consider  a  function 


em  =  f{?)  -  arB  g(?,  5)  f  0<C<1>. 


.  (2.  5.  19a) 


Then 


1^111.  =  0  give8  (ara)min  =  f'(  ?)/g'<  ?,?)  (0  <  5  <  l)  (2.5.  19b* 

where  primes  denote  differentiation.  And 

-LT-  =  f"  (  n  -  (ar8)min  g-  (?,?)  =  g*  (  ?)  [f  (?  )/g*  (  ?,  *  )1 


=  *'««)  l  -3T  '".'min  ’  «  0  <  <  <  *> 


(2.5.  19c) 


The  sign  of  02C/&?2  therefore,  depends  on  the  variation  of  (ar  )  .  as  a 

s  m  l  n 

function  of  $  provided  g'  {?,  5)  i*  shown  to  be  positive  throughout 
0  <  ?  <  1  . 

In  view  of  this  lemma,  a  plot  of  (KC'/tj.)  versus  f  for  variouf 
5  is  given  in  Fig.  4.  It  i3  easy  to  show  that  g'  (  ?»  €  )»  which  corresponds 


fob  l<o  »,  the  value  of 

AT  5«1  li  SMALLER  THAN  AT  'S-O 
AND  TH«  CURVES  A  AC  MONOTOMCAU.V 
DECREASING  WHILE  TO  A  Tf  >0-9 
THEY  ABE  WCRE  ASiNG,  THE  VALUE 
OF  KO'/Etr  AT  5-1  LYING  AEOVE 
THAT  AT  J*»0 


^  -axa 


5-00 

(STONER  LIMIT) 


;  ol  &r  0-4  Oi  r° 

FIG.*:  PLOT  OF  X0/€f  VS.5  THAT  MAKES  ^  )  VANISH  FOR  VARIOUS  A 


\  * 
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a&pa-ii 


to  the  derivative  of  the  exchange  energy  with  respect  to  f  in  0  <-f  <  1 
positive,  by  just  examining  it  directly.  The  plots  in  Tig.  4  show  that 
(K^  monotonically  increasing  with  f  for  {  «s  down  to 

5  =  1,  and  becoming  monotonicall 


y  decreasing  for  ■£  <  0.  9,  the  5  *  0. 9 
curve  being  almoet  flat.  Thie  shows  that  for  0.  9  <  t  <  «o‘,  •  C/lf  >  0  for 
0  <  C  <  1  giving  the  possibility  of  a -minimum  in  this  range,  while  for 
0  <  {  <  0.  9t  <  0  ,  and  hence  the  intermediate  states,  0  <  5  <  1,  are 

relative  maxima  and  so  are  inadmissible  as  ground  states.  After  this,  one 
must  compare. the  interior  minimum  with  the  values  at  tho  end  points  to 
arrive  at  the  absolute  minimum. 


In  the  Coulomb  and  in  the  Stoner  case,  the  following  situations  arise 
and  are  here  given  for  the  sake  of  completeness. 


0:  Coulomb  gas  (Bloch  [3]) 


8  min 


(ar3)m^n  is  a-  monotonically  decreasing  function  of  ?  ,  taking  a  value  of 
about  v  near  ?  =  0,  and  0.  63ff  near  f  =  i.  Hence,  the  intermediate  states 


V  . 
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MWif  V 
-  w  . 


k:._  •  *  * 

■-  ■  '• 
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do  aot  occur  a*  ground  state*  in  the  Coulomb  case.  More  explicitly  from 

(2.  S.  US 

. 

{82£/*f2){or  )  *  -  U*  >  t  1*-y)Z/3s  x2yZ<  0(0  <  £  <  1)  (2.  5.  22) 

•  min  “ 


- 


V . 


T 


,  So  then  on*  need*  to  compare  juat  the  energies  of  the  P  and  the  F  state*. 
Thi*  give*  the  usual  condition*  of  Bioch,n&mely,  fcr  ar#  >  2.839  ,  ona 


ha*  F  . 


To  see  explicitly  the  equivalence  of  (2.  5.  21)  with  consistent 
definitions  of  kpj  and  kpj  f  we  have 

2  ,  ar«  >  or 
M  =  /2m  -  (  -y )  kp  kpj  =  kp.j  /2m  -  (-™*  )  kp  kpj 

from  which  the  rest  are  obvious. 


'  oo  :  Stoner  Gas  (Stoner  [4}  ;  Suri*  [17]  ) 

From  (2.  5:  14),  slnce^{  z  )  ~  l/9  z2  for  z  ~  ao  , 


[(U;i5/3^  (.  ■  ;>5/l  j  --r  (^l(ltf2) 


(0  <  £  <  1)  (2.  •>.  23) 


From  (2.5.16),  since  G(z)~  l/6z^  fo 


Kg’ )  a  n 

CF  min  _ 


(2.  5.  24) 


(K0’/eF)rnin  is  monotonically  increasing  with  £  ,  taking  a  value  of  about 
0.  67  near  £  =  0  and  about  0.  79  near  £  =  1,  and,  hence,  the  intermediate 
states  can  be  ground  states.  More  explicitly  from  (2.  5.  18). 


1 
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(*  -  Y) 

FV  Txy 
F/mln 


(0  <  ?  <  1) 


(2.  5. 


This  is  positive  and  so  the  unsaturated  ferromagnetic  states  can  be  ground 
states.  From  these  one  may  infer  the  usual  Stoner  criteria. 

The  definition  of  chemical  potential  here  is 

^  =  kFf/2m  '  4  2/N<1  +  -  r^N<l-?> 

from  which  the  equivalence  of  (2.  5.  24)  with  consistent  definitions  of 
and  kpj  follow. 

3.  Yukawa  Case 

The  extreme  cases  of  ?  =  0  and  £  =  oo  having  been  reviewed,  the 
cases  with  finite  ?rs  will  now  be  discussed.  First  of  all,  the  possibility 
of  intermediate  states  will  be  discussed  in  view  of  the  lemma  quoted  earlier 
(2.5.19c).  In  Fig.  4,  a  plot  of  <K*'/cr>m£n  versus  ?  for  various  values 
of  £  ranging  from  0.  5  to  oo  is  given.  F rom  the  expressions  (2.  5.  14,  16,  18 
one  notices  that  the  usual  Taylor  expansion  for  small  arguments  when  £  ia 
large,  indicates  the  Stoner-type  behavior.  This  sort  of  argument  gives  the 
inequality  that  for  5  »  2^  (1  ♦  t> 1  /3  one  obtains  the  .hort-ranee  behavior 

1  /3  ^ 

°r  ^  <<  ^  ^  one  obtains  the  long-range  behavior.  From 

Table  3,  we  had  earlier  obtained  an  estimate  which  is  a  factor  -y  of  what 
is  expected  as  above.  More  careful  analysis  shows  that  for  $  ~  0.  9  the 
transition  from  the  "long-range"  to  the  "short-range"  behavior  sets  in,  in 
the  sense  that  the  intermediate  states  tend  to  be  disallowed  for  £  <  0.  9  and 
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allowed  for  £  >  0.  9.  This  is  brought  out  clearly,  when  Fig.  4  is  examined 
in  conjunction  with  the  lemma.  o 

/ 

In  Figs.  5a-e  a  plot  of  C(?)  versus  £  for  various  .£  but  for 

several  fixed  ( KB'  /e F)  are  given.  These  plots  show  some  interesting  features. 

In  Fig.  5a,  KB' -  0.  2  which  in  the  Stoner  limit  exhibits  the  P  state  as 

ground  state.  It  is  seen  that  the  P  state  is  the  absolute  minimum  for  all  5* 

In  Fig.  5b,  K O'  /e^.  =  0.  667,  which  in  the  Stoner  limit  shows  the  absolute 

minimum  precisely  at  ?  =  0.  As  §  is  decreased  the  curves  look  exactly  as 

in  Fig,  5a,  showing  that  for  finite  £  ,  P  ia  still  preferred.  For  KfJ'/e.-,  =  0.  68 

the  Stoner  curve  shows  a  minimum  at  £  =  0.  5,  i.  e.  ,  UF  is  the  ground  state. 

This  minimum  moves  to  £  =  0  quickly  as  £  is  decreased.  This  movement 

of  the  minima  as  £  decreases  is  well  brought  out  in  Figs.  5d,  e.  In  Fig.  5d, 

=  0.  794  where  the  Stoner  limit  gives  the  F  state  a3  the  absolute 

minimum.  On  this  plot,  £  =  6  is  seen  to  exhibit  an  absolute  minimum  for 

some  £  /  1  ,  and  for  £=2x2'  ,  Pis  the  ground  state,  showing  that 

as  £  decreases  the  absolute  minimum  £  moves  to  the  left.  In  the  last 

min 

figure  of  this  series,  K0'/eF  =  2  where  the  Stoner  curve  is  well  in  the  F 
state.  As  £  decreases,  one  gets  UF  for  certain  £  ,  and  for  smaller  £  , 
the  P  state. 

To  show  these  more  explicitly,  a  plot  of  (Ktf'/cp-)  versus  £  is  given 
in  Fig.  2b,  such  that  this  movement  of  the  £  .  *3  brought  out  clearly.  For 

a  given  (KO'/e^),  as  £  is  varied,  we  have  the  minima  at  various  magnetizations. 
Conversely,  for  a  given  £min  , there  is  a  set  of  optimum  values  of  (K0'/«F) 
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which  minimize  the  total  energy.  This  is  obvious  from  Figs.  5a  -  e. 

♦  co,  above  K0'/*F  =  0.  79  one  has  F  state  and  below  K Q'/t-v  =  0.  67 
one  has  P  state  and  between  these  two  one  has  UF  as  ground  states.  As  £ 
is  decreased,  the  region  of  UF  becomes  smaller  and  smaller,  giving  a  higher 
and  higher  (K0'/eF)  and  finally  at  about  £  =  0.  9,  coalesce,  showing  that  UF 
is  no  longer  preferred  below  £  =  0.  9.  Below  this,  one  uses  the  condition 
obtained  by  just  comparing  the  F  and  P  states.  For  the  sake  of  showing 
how  the  Bloch  limit  is  obtained  for  £  =  0,  a  similar  plot  of  fn  (ar()  versus 
5  is  given  in  Fig.  2a.  On  this  same  figure,  a  plot  of  the  Thomao-Fermi 
screening  versus  In  (arg)  is  also  given.  This  plot  never  crosses  the  ferro¬ 
magnetic  region  for  any  (arg),  showing  that  F  state  never  occurs  in  HF 
for  an  electron  gas  with  Thomas-Fermi  screening. 

(b)  Just  for  the  sake  of  completeness  we  derive  the  expression  for 
<H>/CIq  in  the  Overhauser  SDW  case.  From  (2.  5.1),  using  the  Overhauser 
Green  s  functions  (2.  3.  1)  and  (2.  3.  6)  or  evaluating  H  ^  directly,  we  get 


For  £ 
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<H>  f  d3q  [7 

^  li~ 


a2  2*  <a+9>2  2 

[lisr coa  V  ~~*m  8in  *a 


)  nF<u+^a»  + 


/.a2  2  <a+9>2  > 

+  [  *a  8 in  0  +  — -a": cos  6 


,3  ,3 


2m 


)  °F  (u.(a»l 


(cos  0  nF(u+(a))  + 


.  2 


+  Bin  0  np-(u_(a)))(cos  0  , nf.(y+(a'))+ ,nF(*.»  (q')))  + 


.  2 


+  (3in^nF(ij+(a))  +  co3^nF((j_(a)))  (sin^.n^u^'JJ  +  cos^.n^u (a')))  + 


/gt  i  (a)git(a') 

+  2/g]'ifa^n<a'>  cos 0asin0aco^a,3in0a,(nF((J+(q))-nF(tJ_(a)) 

(nF(u+(a'))-nF(w  (a')))l  (2.  5.  26) 

If,  as  Overhauser  [18]  did,  only  the  u+(q)  branch  is  assumed  to  be  occupied 
and  further  gj  J  =  gj  j  =  g,  we  recover  his  expression.  Incidentally,  this 
generalizes  Overhauser's  results  to  finite  temperatures.  Treating  0^  as 
a  variational  parameter,  one  recovers  back  the  condition  on  tan  20  which 

a. 

can  be  obtained  from  (2.  3.  8)  also.  Moreover,  this  generalizes  Overhauser's 
results  for  intermediate  magnetizations  whereas  Overhauser's  work  concerned 


only  the  paramagnetic  state. 

Summarizing,  we  have  in  this  chapter  discussed-single  particle 
states  and  the  ground  states  of  a  polarized  electron  gas.  The  discussion 
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includes  some  formal  results  about  SDW  (symmetry  breaking  solutions) 
and  electrons  in  a  periodic  potential  (Bloch  electrons).  It  is  shown  that  the 
exchange  energy  in  general  is  a  very  complicated  function  of  ?2  ,  and  not 
just  a  power  series  in  ?2  .  The  main  results  are  displayed  in  Figs.  2a,  b 
which  give  the  regions  of  magnetizations  for  optimum  values  of  the  coupling 
strengths  and  screening  which  minimize  the  total  energy  of  the  system.  The 
use  of  the  Thomas-Fermi  screening  is  valid  under  certain  approximations 
for  a  Coulomb  gas,  and  implies  a  definite  relationship  between  the  screening 
and  (arg)  .  From  this  it  is  shown  that  in  HF  the  electron  gas  can  never 
become  ferromagnetic.  All  the  calculations  reported  here  are  within  the  HF 
approximation  and  for  zero  temperature,  though  the  formalism  contains 
results  for  finite  temperatures  also. 


III.  INDIVIDUAL.  PARTICLE  AND  COLLECTIVE  EXCITATION  OF  THE  SYSTEM 
1.  General  Introduction 

In  order  to  study  the  collective  excitations  of  the  system  one  needs  to 
study  the  characteristic  correlation  functions.  The  formalism  developed  in  the 
last  chapter  can  be  used  for  this  purpose.  In  this  section  the  necessary  equations 

i 

are  derived  for  the  correlation  functions  of  interest.  In  view  of  the  complications 
drastic  approximations  are  made  to  solve  them.  These  will  be  discussed  in  the 
next  two  sections. 

The  first-order  response  of  a  system  to  a  given  external  perturbation 
can  be  expressed  in  terms  of  retarded  correlation  functions.  There  is  a  locus 
of  poles  in  these  functions  in  the  frequency-momentum  plane.  This  locus  gives 
a  relation  between  the  frequency  and  the  wave  vector  corresponding  to  the 
forcing  field.  These  are  the  bound  states  or  the  collective  states  of  the  medium 
indicating  resonant  response  to  the  external  field.  There  is  a  second  frequency- 
wave  vector  relationship  which  is  not  a  pole  structure  in  the  correlation  function 
but  which  is  a  branch  singularity.  The  branch  singularities  correspond  to  the 
scattering  states  of  the  system.  In  general,  the  poles  of  the  correlation 
function  have  a  real  and  an  imaginary  part.  The  locus  in  the  frequency-wave 
vector  plane  which  makes  this  imaginary  part  non  zero  gives  a  locus  of  a 
complex  pole  in  this  plane,  showing  that  the  collective  part  has  a  finite  life¬ 
time.  When  this  imaginary  part  is  zero  one  has  a  real  pole  which  is  the  bound 
state  excitation  which  occurs  only  for  a  certain  range  of  the  wave  vector  in  a 
certain  approximation.  These  may  be  stated  physically  as  follows.  In  looking 
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for  the  collective  excitations  normally  examined  (and  also  those  studied  here), 
one  looks  for  coherent  motion  of  a  particle-hole  pair  with  certain  characteristics, 
like  their  spin,  wave  vector  and  energy.  The  bound  state  corresponds  to  this 
coherent  motion.  Beyond  a  certain  value  of  the  wave  vector  (difference  in  the 
wave  vectors  of  the  pair)  which  is  the  wave  vector  of  the  collective  mode,  this 
coherence  is  lost  and  one  no  longer  has  a  collective  mode.  This  is  the  beginning 
of  the  scattering  states  of  the  system.  Thus,  the  nature  of  the  singularities  in 
the  retarded  correlation  functions  tells  us  about  the  nature  of  the  excitations 
in  the  medium.  Here,  in  particular,  the  responses  of  the  system  to  three 
kinds  of  external  fields  are  studied:  (1)  oscillating  fields  giving  rise  to  the 
fluctuations  in  the  charge  density  of  the  medium  thus  yielding  the  density-density 
oscillations;  these  may  also  bethought  of  as  a  singlet  state  oscillation  of  the 
particle-hole  pair;  (2)  oscillations  of  the  z  component  of  spin  density  by  means 
of  an  external  magnetic  field  in  the  z  direction  fluctuating  the  z  component  of 
spin  density  so  as  to  create  spin  density  oscillations  without  the  attendant  flipping 
of  spin;  these  are  the  triplet  oscillations  of  the  pair  with  their  projection  zero; 
and  (3)  fields  which  fluctuate  the  transverse  components  of  spin  density  by  means  of 
an  external  magnetic  field  in  the  othe r  transve rse  direction  giving  rise  to  spin 
flip  oscillations  but  without  the  fluctuations  in  the  total  charge  density,  thus 
leading  to  spin-wave  type  collective  modes.  These  are  the  triplet  state 
oscillations  of  the  pair  whose  projection  is  unity.  Instead  of  examining  the 
retarded  correlation  functions  themselves  directly,  a  set  of  time-ordered 
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functions  which  are  related  to  these  (in  the  same  way  as  in  Chapter  2  concerning 
the  single  particlT'sTatc*.}_will  be  studied.  This  is  done  because  the  variational 
derivative  technique  can  be  used  to  advantage. 

Define 

•®|(  1 )  =  <rt(  1  )  -  <<r.  >  (i  =  0,  1,  2,3)  (3.1.1) 

where  <r.  is  the  operator  tr  ( (// +  7^  ^ ) ,  \}/  is  taken  to  be  a  1x2  column 

matrix,  and  <  o\  >  is  the  usual  average  density  as  in  (2.  1.  21a,  b),  where 
stands  for  the  Pauli  matrices  and  r  for  unit  matrix.  Thus,  <r  (  1  )  =  excess 
particle  density  over  the  mean  particle  density,  p(  1  )  ,  £  {  1)  =  excess  spin 
density  over  its  mean.  With  these  definitions  the  following  correlation 
functions  are  computed: 

X-jOl')  =  1  )  a.  (  1 1 )  )  > 

(»,  J‘  s  0,  3)  •  (3.  1.2) 

The  time  ordering  here  corresponds  to  the  Bose  operators  in  contrast  to 
Chapter  II.  For  i  =  j  =  0,  (X  QQ  =  )  this  function  is  related  to  the  dielectric 

constant,  S  ;  for  i  =  j  =  3,  it  is  related  to  the  longitudinal  spin  susceptibility, 

^  33  ’  ^or  *  :  M  +  ),  j  :  2  (  -  )  ,  it  is  related  to  the  transverse  susceptibility 
X+  .  (1,2  now  refer  to  the  spherical  components  (  ( cr j  +  i  <r^)/2)  instead  of 

Cartesian  components  as  in  Chaper  II  ). 

As  in  Chapter  II,  Section  1,  a  set  of  formal  relationships  can  be 
derived  just  on  the  basis  of  the  definitions  (3.  1.  2).  These  are  given  in 
Appendix  A  . 
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To  relate  the  functions  (3.  1.2)  to  the  variational  derivatives  with  respect 
to  external  fields,  let  us  first  note  that  the  fields  U|  and  Uj  introduced  in  the 
last  chapter  may  be  combined  to  give  the  generators  of  total  density  and  the  z 
component  of  spin  density  by  forming  the  following  combinations: 

Uj  =  Uq  +  U3  and  Uj  =  Uq  -  Uj  (3.  1.4) 

so  that  (2.  1. 22)  may  now  be  written  as 

H'  -  j*  d3i  [UQ{it)  n ( rt)  f  U3(rt)  <r3  (rt)]  (3.  1.  5) 

where 

n  (rt)  =  i//j+  (rt)  if/  j  (rt)  +  \f/^  (rt)  l//j  (rt) 

and 

<r3(rt)  =  xfy  j"  (rt)  l//|(rt)  -  j"  (rt)  l//j  (rt) 

In  view  of  this  we  have 

6,5  6  6  6  6 

"STJj"  7>  Uj  STT* »  Tuj~  '  (TUj  “ 

To  generate  o-+(rt)=i//^  (rt)  i//  j  (rt) 
and 

cr_(rt)  =  (rt)  i//|  (rt) 

one  introduces  the  following  additional  Hamiltonian 


H"  =  ^  d^r  [U+  (rt)  a  (rt)  +  U  (rt)  a+  (rt)] 


(3.  1.6) 


The  fields  (U+  U  U3),  when  suitably  chosen,  are  the  components  of  a  space- 
and  time-dependent  external  magnetic  field.  The  equation  for  the  single  particle 
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Green's  function  in  the  presence  of  all  these  fields  may  now  be  written  down  as 
in  the  last  chapter  (here  no  new  notation  is  introduced  to  distinguish  this  from 
the  old  one  as  this  does  not  cause  any  confusion)  corresponding  to  (2.2.  28) 

V  ^ 

G‘1(11,)  =  [iiq  +  TST  "  v(D-Uo(D-  r_u+(i)  -  r+uji) 

-  t.3  u3(i))  6(4,.(1-1') 

+  i  j*  d4T2/(l-T)  [tr  G(TT+)]  6{4)  d-1') 

tijd4!  d4r^(i-DG(iT')r  (Pi'jD  .  (3.1.7) 

To  make  use  of  the  variational  derivative  technique,  the  expressions  (3.  1. 2) 
will  now  be  recast  in  terms  of  suitable  derivatives.  For, 

<<r.  (1)  >  =  -  i  [tr  t^G  ( 1 1  *)]  (i  =  0,  1,  2,  3) 

where  1,  2,  3  now  refer  to  spherical  components  (+,  3).  So  one  has, 

-  -  gu^T-TT-  [tr  TiG(ll+)]=  <T(o\(l)oj(2))>  -  <  <r.(l)>  <<r.(2)  > 

=  <T((a.(l)-<a.(l)>)  (a.(2)  -  <a.(2)>))> 

=  <T(<r.(l)  ^(2))> 

Al30 


6  G  ( 1  1  ) 


^  d42  d^GV-Z)  T.  (2  3;  2)  G(Tl  (tj+)  ) 


6  U.  (2) 

J 


(3.  1.9) 


t 


k  ' 
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from  (2.  1.  25),  and^defining 


-1, 


r.(  12 ;  3)  =  6G~‘(12)/6U.(3) 

J  J 


(3.  1.  10) 


The  r.  are  matrices  like  T  in  the  last  chapter.  Hence 

i  ° 


Xj.dZ)  =  i  [  d4Td42  (tr^GdDr.dZ^)  GdUtj1)  )}]  .  (3.1.11) 


The  equations  for  I\  can  now  be  written  down  with  the  help  of  (3.  1 . 7) 


r (  1  2 ;  3 )  =  -Tj  6  (4)  (1 -2)  6 (4)(1  -3)  - 


-  i 


j*  d4Td4Id45^(l-T)  [tr  {G(ri)r.(2l;  3)  G  (  5T  ( t/)  )} ]  6  (4)(1  -2) - 


-  i  j*  d4T  d4l'  d42  d432/(l-l)  G(1  2)r.  (15;  3)  G  ( IT')  Tq  (  T' 2;  T)  + 

fir  (T2;  n 


r  a_  a  _  um1 

+  i  \  d4l  d T'  2/(l-l)  G  (11*)  -gq- 

J  j 


nr 


(3.  1.  12) 


These  equations  are  exact.  From  (3.  1.12)  one  can  derive  the  equations  for 

6ro/6  Uj  and  so  on.  These  form  a  set  of  linear  equations  for  Tj  ,  F^  ,  Fj 

except  that  these  are  all  coupled  to  T  which  satisfies  a  nonlinear  equation. 

*  .  "> 

The  second  term  in  (3.  1.12)  has  an  obvious  interpretation  and  lets  us  introduce 


7.(12)  =  i  j'd4Vd4Z[tr[G(lT)r.(TI;  2)  G(21  (t^))]  ]  (3.1.13) 

This  term  arises  from  the  Hartree  term  in  G  ' ,  the  next  to  last  term  in 
(3.  1.  7),  and  has,  therefore,  the  interpretation  that  it  is  the  response  of  the 
total  particle  density  to  the  component  IT  of  the  external  field 


i 


1 
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yxlli)  *  -f <T  <p(D, ^(2)  )>  .  (3.1.14) 

,  The  various  terms  in  (3.  1.12)  may  be  interpreted  as  follows.  The 
first  term  is  the  lowest  order  HF  term  since  in  the  last  chapter  this  term  was 
.seen  to  give  the  conventional  HF  single  particle  excitations.  The  second  term 
is  due  to  direct  scattering  processes  involving  fluctuations  in  the  Hartree  self 
energy;  the  third  term  arises  from  the  exchange  scattering  processes  and  the 
last  term  includes  all  the  higher-order  vertex  corrections  coming  from  various 
other  scattering  processes.  The  Eq.  (3.  1.  12)  in  its  full  generality  is  very 
complicated  and  simplifying  approximations  have  to  be  made.  These  will  be 
discussed  in  the  next  section  as  well  as  the  results  obtained  therefrom. 

2.  Approximations  (Spatially  Uniform  Solutions) 

We  now  consider  the  free  gas  where  the  single  particle  potential  is 
taken  to  be  zero.  Then  we  have  a  spatially  uniform  system  where  the  Green's 
function  G(ll')  depends  only  on  the  difference  ( 1  - 1  * )  and  the  functions 
I1.  (  1  2  ;  3  )  depend  only  on  (1-2)  and  (1  -3).  Before  proceeding  further,  we  note 
that  the  last  term  in  (3.1.12)  is  of  order  as  can  be  seen  by  forming  the 

equation  for  .  We,  therefore,  drop  this  term  in  all  our  calculations. 

This  amounts  to  neglecting  all  vertex  corrections.  (See,  for  example, 

Rajagopal  [45]  for  a  discussion  of  these  equations  in  the  extended  RPA  scheme 
using  the  above  formalism.)  We  now  take  the  Fourier  transform  of  this  using 
(2,  2.  1)  and  (2.  5.  4)  where  4-D  notation  is  used  for  convenience. 
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r.(k;  q)  =  -  t.  -^UW^q)  -  i  ^q^p  ( q+ q)  •• 


r.  (q;  q)  G(q)  T  (k;  q-k)  -2^- 

1  (2») 

y.  (q)  is  just  the  Fourier  transform  of  (3.  1.  14)  and  has  the  form 

y.  (q)  --  i  f  d  [tr  f  G(q  + q)  T  (q;  q)  G(q)l  1 
;i  J  (25)* 


(3.  2.  1) 


(3.  2.  2) 


From  the  structure  of  (3.  2,  1),  following  a  suggestion  used  in  another  context 
by  Nambu  [43]  ,  an  algebraic  transformation  may  be  made: 

F(k;  q)  =  At(k;  q)  +  ^(q)  Ti(q)AQ(k;  q)  (3.  2.  3) 

where  now 

\(k;  q )  *  -T.  - 


-  i  {  2/(  k  -  q  )  G  (  q  +  q  )  A.  (  q  ;  q  )  G  (  q  )  To  (  k  ;  q  -  k  ) 

J  [my 

In  view  of  (3.  2.  3)  ,  one  has  the  form 

/-  4  — 

7.  (  q )  =  i  J  d  ^  [tr  |  G  (  q  +  q  )  Ai  (  q  J  q)  G(q)^  ]/C(q) 


7  (  q )  =  i  i 

*  ,4— 

1  d 

1  (2*)4 

E(q)  =  1  - 

-i^(q)  | 

d4  q 
(2jt)4 

[tr  |  G  {  q  +  q  )  Aq  (  q  i  q)  G(q)^  ] 


(3.  2.4) 


(3.  2.  5) 


(3.  2.  6) 
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But  7Q(q)  is  the  Fourier  transform  of  the  density  correlation  function,  and 
so  one  has  the  following  general  relationship  between  E(q)  and  7Q(q)  which 


is  obtained  from  their  definitions  (3.  2.  5,6)  : 

Vq)  1  37T W  [  TrTTT  '  1  1 


(3.2.7) 


Comparing  with  (2.  50)  of  Pines  [32]  (page  41),  one  thus  sees  that  E  (q)  may 
be  identified  with  a  propagating  dielectric  constant  of  the  medium.  In  view  of 
the  definition  (3.  2.  6),  the  equation  for  T  (  cf.  (3.  2.  3))  may  be  rewritten: 


(3.  2.8) 


To(k;  q)  =  Ao(k;  q)/E(q) 
and  the  equation  (3.  2.  4)  for  Aj  : 

A  j  (  k  ;  q  )  =  -  Tj  - 


~  i  j*  r/4k-:k)~  G(q  +  q)  A.(q;  q>  G(q)  Ao(k  ;  q-k)  .(3.2.9) 

Then  one  finally  arrives  at 

X33  =  ‘  y  ^  d  ^  ftr[T3  G(9+9")  A  jtq-;  q)  G(q)}  ]  + 

+  «  ^(q)  73(9^  d  ^4  [tr  |r3  G(q+q)  Aq  (q;  q)  G(q)|  ]  (3.2.10) 

x  +  -  =  *  5"  ^tr  (Ti  G^9+q)  A 2  ( q ;  q )  G  ( q)^  1  + 


+  i  V(  q)  7  2  (  q  )  ^  ( tr  {  Tj  G  (  q  +  q  )  A  Q  (q  ;  q  )  G  (  q  )^  ]  .  (3.  2.  1 1) 
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In  this  technique  the  dielectric  screening  of  the  exchange  term,  namely  in  the 
second  term  of  (3.  2.  9)^  cornea  about  quite  naturally  without  any  ad  hoc 
assumptions.  Moreover,  a  perturbative  type  of  solution  of  the  Aj  equation  in 
powers  of  the  interaction  strength  gives  the  results  of  the  extended  RPA  [45]  . 
We  will  now  present  four  cases,  three  of  which  depend  on  further  simplifying 

approximations  of  the  equations  {3.  2.  9). 

1 

To  motivate  the  approximations,  let  us  examine  how  the  terms  in 
(3.  2.  9)  came  about.  To  this  end,  the  expression  (2.  5.  5)  for  G  *  (k)  with 
all  the  =  0  ,  must  be  seen.  This  contains  a  term  involving  r  ,  or 
equivalently  Aq  .  It  is  this  which  couples  the  various  R  to  Rq  as  well  as 
making  the  equation  for  R^  nonlinear.  However,  the  Rq  equation  Is 
inhomogeneous.  It  was  earlier  seen  that  taking  the  inhomogeneous  term  for 
rQ  to  be  the  solution  gives  the  HF  approximation  for  G  .  So,  if  the  equations 
for  R.  or  equivalently  for  A^  are  derived  from  G^^.  ,  one  will  arrive  at 
equations  of  the  form  (3.  2.  1),  but  in  the  last  term,  Rq  appearing  at  the  end 
is  set  equal  to  -1  thus 

r.  (k;  q)  =  -  rt  -V(q)  7^  (  q  )  + 


i  T  -^-3“.'2/{k- q)  G{q  +  q)  r  (q;  q)  G(q)  . 
J  ^  1 


+  i 

w 

And  correspondingly 

Aj  (k;  q)  =  -  t.  + 

,4- 


(3.  2.  12) 


+  i  f S'(k-q)  G(q  +  q)  A  .  ( q  J  q)  G(q)  , 

J  (27rr  1 


(3.  2.  13) 
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These  equations  then  constitute  the  usual  RPA  equations  but  include  lowest 
order  exchange  terms  and  must  be  solved  to  arrive  at  the  various  response 
functions.  Again,  the  equation  for  A.  as  in  (3.  2.  13)  cannot  be  solved  completely, 
but  this  time  for  a  different  reason.  The  lowest  order  RPA  consists  in  taking  the 
inhomogeneous  term  to  be  the  solution  for  .  But  this  is  the  lowest  order 
approximation  when  k)  depends  on  k  .  On  the  other  extreme,  if  ^/(k) 
is  chosen  to  be  independent  of  k  ,  which  is  true  for  zero  range  interactions, 
then  further, if  G's  are  assumed  to  be  diagonal,  (3.  2.  3)  can  be  solved  trivially 
for  .  This  then  is  the  second  case  where  some  physically  nontrivial 
results  are  obtained.  In  this  way  of  treating  the  problem,  it  is  thus  seen  that 
the  results  for  the  Coulomb  gas  and  the  zero  range  gas  come  about  quite 
trivially.  These  will  be  discussed  in  the  next  two  subsections.  We  henceforth 
assume  the  G's  to  be  diagonal. 

There  is  another  possibility  in  the  present  context.  This  is  to  consider 

the  full  set  of  equations  for  A.  and  G  *  together.  We  assume  A  =  -1 

i  o 

to  start  with.  This  is  substituted  in  the  expression  for  C(q)  which  now 
becomes  the  RPA  dielectric  constant.  This,  therefore,  screens  the  exchange 
term  in  G  *  .  The  resulting  equations  for  A  .  ( i  =  1  ,  2  ,  3  )  are  now  independent 
of  Aq  but  are  of  the  form  (3.  2.  13)  where  in  the  last  term,  the  interaction 
potential  is  screened  by  .  This  new  set  will  be  discussed  to  justify 

the  use  of  the  Yukawa  potential  in  this  problem.  This  is  discussed  as  Case  (c). 

After  describing  these  results,  in  the  next  section,  the  results  are  extended 
to  Bloch  electrons,  which  form  the  fourth  case  of  the  present  discussion. 
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It  may  be  mentioned  in  passing  that  the  same  method  was  employed  by  the  author 
[46]  to  derive  quite  trivially  the  dielectric  constant  of  an  electron  gas  in  the 
presence  of  an  external  magnetic  field. 

A.  Coulomb  Gas 

f  4ir^z 

Here  (q)  =  — —  and  only  the  most  trivial  approximations  can 

q 

be  made  to  give  closed  form  results.  Here  the  solutions  to  (3.  2.  13)  are 
taken  to  be 


q)  =  -  Tj 

Then 

Sjtq)  =  1  +  V  (q)  (  Aj  +  Bj ) 

i  r 

To  (C1)  =  *  B-y 

Y^q)  =  o 

72r(q)  -  0 

73!(q)  =  - 


X+r  (q)  =  -  Cj 


(Ar  -  Bj) 

1  +tMq)  (A  j  +  Bj) 

Aj  +  Bj  +  4<y(q)  Pl  Bj 
i  +  °i$  (  q )  (Aj  +  Bj) 


(3.2.  14) 


(3.2.  15a) 


(3.  2.  15b) 


(3.  2.  15c) 


(3.  2.  1  5d) 


(3.  2.  1  5e) 


(3.  2.  1  5f) 


(3.  2.  1  5g) 
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Here 

■  I  r> 

AI  =  1  J  G||  (l  +  «)G|t  «*> 
Bi  *  >1  GH  <’+’>Gtl<’> 

ci  =  ‘I  f^fpr  Gtt  <<)^>Gti  <«> 


(3.2.  16a) 

(3.  2.  16b) 


(3.2.  16c) 


Integrals  of  the  type  (3.  2.  16a,  b,  c)  with  the  same  general  structure  will  appear  in  the 
future  also.  Using  a  method  given  in  KB,  these  can  be  evaluated  in  general  in 
the  form: 


A  (  £{J  )  =  ^ 
B { )  =  J 
C  { <ju )  r  y 


d3a  i 

r 

r  dtj 

(2  7T)3  ' 

'  2  IT 

J  2» 

d3a 

f  du 

f  dj 

(27T)3  ' 

'  L* 

^  2jt 

d3a 

f  du 

r  du 

(2^)i  « 

'  2  IT 

^  zv 

Af I  (a+a;  u>Af f  (a;  «*>> 

u  -  U  +  <J 

_  — .  —  S3 

A 1 1  ( q  +  q ;  u)  Aj  |  ( qj  «•>) 


(n^-(u)-n.j.(u)) 


(3.  2.  17a) 


u  -  u  +  U 

Aj  |  { a+  q ;  •*>)  A|  j  ( a; 

U)  -  u  +  w 


(n^i^-n^  (u)) 

(3.  2.  17b) 


(nf  (u)-nF(cj)) 


(3.  2.  17c) 

The  weight  functions  A  ,  are  as  given  in  Chapter  II.  Hie  HF  approximation 
gives 


Aj  |  (  cp)  =  Zu  5  {u  -  e  jHF(  a>  J  Aj  |  {  c^id)  =  2  n6  (  u  -  e|IF  (  q))  .  (3.  2  18) 


These  in  (3.  2.  17)  give  the  combination 

[nF(€<TF  {a+a»  -  nF(e(r^F{a^l 
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of  Fermi  functions.  We  take  this  to  be  ■ 

K(kF«r"  la  +  al>  -  n+(kF<r,  -  III)] 

for  T  =  0°K  as  in  (2.  5.  8).  Note  that  these  integrals  in  general  have  a  principal 
part  and  a  delta  function  part  or  in  other  words  they  have  real  and  Imaginary 


parts.  For  the  present  case: 


^  =  f  — 

1  J  (2JT) 


r  d  a 

Bl  =  \  - T 

1  J  UTt) 

.3- 

r  d  a 

ci  =  \  — r 

1  J  (2tr)J 

u+!  (  q)  =  qZ/2m  - 


nF  *  a + a  ®  -  nF  <a» 

u  -  w+  ( a  i  5)  +  w+  (a) 


( a+ 5))-nF  a) ) 

- 1 - z - rz — 

u  -  u_  ( a+ a) +  u.  ( a) 

nF(u_r(a  +  a))  -  nF(w+I(g[)) 

- T - z - T“z - 

U  -  u_l( a  +  a>+  ( a) 


-(4) 


qkFa  +  2 


1  ,  2  .  2  ,  q'kFa 

J  (q  ’  kFa  >  ,n  1  q  +  k  ■■ 


Here  kF^  :  +  corresponds  to  kf|  J  "  to  M 


(3.  2.  19a) 


(3.  2.  19b) 


(3.  2.  19c) 


(3.  2.  20) 


The  expression  for  £j(q)  derived  under  these  approximations 
corresponds  to  the  dielectric  constant  derived  by  various  authors  within 
the  RPA  (see  Pines  [32]  )  in  the  unpolarized  limit.  Wc  will,  therefore,  call 


this  plasma  RPA. 


/ 
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One  notices  from  the  form  of  the  Eq.  (3.  2.  13). and  the  solution  assumed 
(3.  2.  14)  that  the  exchange  terms  are  entirely  neglected  in  this  plasma  RPA. 
Because  of  this  the  polarized  system  behaves  as  though  one  had  just  two  non¬ 
interacting  plasmas  and  the  medium  shows  no  spin  wave  mode  in  X ,  . 

*T  *• 

However,  on  taking  the  exchange  terms  into  account  by  examining  the  T, 
equation  in  more  detail,  as  will  be  done  later,  spin  waves  do  result.  (As  a 
check  on  the  calculations,  if  the  system  were  unpolarized,  A  =  B  =  C,  the  above 
results  reduce  to  the  known  results  found  in  the  literature.)  In  this  approximation, 
Yl(<l)  =  'y2^t0  =  ®»  th*3  implies  that  there  are  no  density  fluctuations  caused 
by  the  fluctuations  of  the  transverse  components  of  the  internal  polarization. 
However,  /  0  ,  showing  that  the  longitudinal  component  does  affect  the 
density.  These  results  are  reasonable  on  physical  grounds  also.  In  the 
complete  F  state,  B  =  Q,  (J  state  is  empty)  there  is  still  a  plasma  type 
excitation  found  in  both  the  density  response  and  in  X33.  (The  complete  F  state 
corresponds  to  an  electron  gas  with  one  spin  state  occupied;  two  particles 
cannot  come  to  the  same  spot  by  dint  of  the  Pauli  principle.  This  is 

significant  as  will  be  seen  in  case  (B).) 

As  a  further  crude  approximation,  to  illustrate  the  results  obtained  so 
far,  if  the  exchange  terms  in  (3.  2.  20)  are  expanded  for  small  q  ,  one  obtains 
an  effective  quadratic  dispersion  for  the  single  particle  energies  with  some 
effective  masses: 


u+(^)  =  q2/2m  - 


2eZk 


F  tr 


with 


1 


m 

IT 


n 


(3.  2.  21) 
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This  is  only  to  give  a  heuristic  argument.  A  more  useful  but  perhaps  irrelevant 
approximation  is.  just  to  assume  a  quadratic  dispersion  law  for  the  two  states  with 
effective  masses,  but  without  the.  significance  of  (3.  2.  21).  In  any  case,  one  has  ; 
plasma  oscillations  in  the  medium,  appearing  in  the  density  response,  yQ  , 
Yj(q)  ,  and  in  .  In  the  P  state,  the  plasma  mode  appears  only  in 

Vq)- 

We  now  derive  the  Thomas-Fermi  screening  in  the  polarized  medium, 
using  an  effective  mass  approximation,  for  the  f  and  l  bands.  F  or  this  we 
need  to  study  C(q,  u  =  0)atT  =  0°K  for  q  ~  0  .  We  then  have 


G(  q.  w  =  0) = 1  + 


n+(kF|  -  la+al>- n+(kF|  -lal) 
I  a+I  lz  * 


+  2m 


I 


d3q  *4<kFj  -li+ali-n+ovi  -  laf> 


(2TT)7  |  q+q|2  -  qZ 


Evaluating  the  integrals  for  small  q  ,  one  has 

2  f  mi  mj  k^-j 


G ( q_,  u=  0  )  =  1  +  ^4- 


‘t  KFf 


2t 


-  1  +  ?ZTF(0  kFZ/q2 


where 


?tF(H 


—  (ar  ) 
V  s 


1  /3 

— L  { l  +  ?)  '  + 


nr 


ml  1/3 

(l-n  ' 

m 


(3.  2.  22) 
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Here  m  are  taken  to  be  arbitrary.  Using  mi  =  mi  *  m,  we  have  plotted  in 

( 

Fig.  2a  1)  versus  In  (ar#)  to  show  where  the  Thomas-Fermi  screening 

lies  on  the  phase  plot  of  HF  ground  states  for  the  Yukawa,  potential.  This  shows 
that  ferromagnetism  does  not  occur  in  HF  for  Thomas-Fermi  screening. 


B.  Zero-Range  Gas  (Stoner  model) 

Here  2/ (q)  =  2/  =  a  constant.  In  this  case,  the  equations  for  Aj  can 
all  be  solved  trivially  if  the  G's  are  diagonal.  These  are: 


A  “(kj  q)  =  - 


1/(1-  2/An) 


0  \ 
1/(1  -  t/Bn) 


(3.  ?.  23a) 


II 


A  j  (k;  q)  =  -  /(l/O-^Djj)) 
A”  ( k  ;  q)  =  -  T*  (1/(1-  2/Cn)  ) 

I  1/(1 -27ad) 


ii 


A/  (k;  q)  -  - 


\ 


C 

whc  re  =  i  J  Gj  j  (  q  +  q  )  Gj  j  (  q  ) 


Then 


II,  .  1 

7o  q  *  ' 


l/d-^Bj) 


/ 


AII  +  BII  "  ZtfAU  BII 


(1  -VAn)  (l-VBn) 


(3.  2.  23b) 
(3.  2.  23c) 

(3.  2.  23d) 


(3.  2.  24a) 
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Xfn(q)  =  -  C„/(l  .7cn)  . 


(3.  2.  24b) 

(3.  2.  24a') 

(3.  2.  24c) 
(3.  2.  24d) 

(3.  2.  24e) 

(3.  2.  24f) 

(3.  2.  24g) 


u"  =  q2/2m  -  n27(  1  +  ?) 

=  q  /2m  -  K6'  (  1  +  ?  )  in  Stoner's  notation  (3.  2.  25) 

As  with  the  results  of  case  (A)  ,  -y^11  ,  -y^11  ,  and  X3  “  show  the  same 
pole  structure.  However,  X.j._  contains  a  pole  also,  showing  the  possibility 
of  spin-wave  modes  in  the  medium  without  additional  approximations.  In  the 
paramagnetic  limit,  An  =  Bn  =  Cn  ;  yj1  contains  a  pole  while  X  ^  has 
none,  for  repulsive  potential.  In  this  limit,  the  results  correspond  to  those 


<v ■  * -at 


t 

L 


■  f 


Ktv  1 
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of  Gottfried  and  Pieman  [30]  and  Wolff  [31]  ,  respectively.  If  one  puts?j>< q)  =*$ 
in  the  results  of  case  (A),  in  the  P  state,  one  finds  it  necessary  to  take 

2/1  q)  ~  to  £ind  “•'A**  ^  agreement  with  those  of  the  present  section. 

This  factor  of  appears  because  the  exchange  contributions  were  neglected 
in  case  (A)  and  were  taken  into  account  here.  For  the  same  reason,  the 
substitution  2/(q)  =  2/  in  case  {A)  does  not  give  the  same  results  as  case  (B) 
in  the  polarized  case.  The  appearance  of  spin  waves  as  a  pole  in  is 

entirely  due  to  exchange.  The  above  results  have  also  been  derived  by  Kubo 
et.al.  [26]  and  Baym  [27].  As  in  (A)  ,  7j  =  T2  =  0  for  the  same  reason. 

The  integrals  can  all  be  formally  reduced  to  quadratures  for  finite  T, 
and  to  closed  expressions  for  T  =  0°K  ,  thus  (only  principal  parts  are  here 


given): 


AII  = 


‘7J?~  J  d<l'  q'  nF<u  n(q'))  /n 
4ff  q  **  c  T 


B 


m 


II 


CII  = 


2  §  d<l'  CII  In 


4j r  q 


^dq'  cl'  j  np-(u_II(q')  fn 


t/_-  (qZ/Zm  -  3^_/rrx}Z 
u  -(q  /Zmfqq'/m) 

/2m  -q  q'/m)2 
£  2  ~  - j- 

u  -(q  /2mfqq'/m) 

+  q2/Zm-q  q'/m 
q  /2m  fq  q'/m 


\ 


II 


-  nF(u+(q'))  In 


ft-  q^/Zm-qq1  /m 


&  -  q  /2mf  qq'/m 


J 


=  (u  -  2K0*£). 


(3.  2.  26) 


/ 
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For  T  =  0°K  ,  these  integrals  can  be  carried  out  and  are: 

| 

m* 


‘II 


3  f 

m  } 

l 


2q3  kF|  /m2  + 
2  2,  2  /  2 


+  0U  +  <1  "  q  /m  ]  in 


(j-t-q  /2m-t-qkp|  /m 


U+q  /2m  -  qkj  ,|  /m  J 


[((J  -  q2/2m)2  -  q2kFj /m2]  in 


B 


II 


'll 


"w  -  q  /2m  +  q  kp|  /m 
W  -  q2 /2m  -  qkpj  /m 

(same  as  above  with  k^.|  replaced  by  k^j  ) 


} 


3 

m 

2~y 


((fi-  q2/2m)kF|  -  (fi+ q2/2m)kf,|  ]  + 

q2/2m+  qkpj  /m 
J^f2+  q2/2m  -  qkf|  /m  J 

iT-q  /2m  +  qkj-j/m 


8*  q 


+  [(^+  q2/2m)2  -  q"  k2  j  /m2]  in 


-r  ((fi-  q2/2m)Z-  q2k2  j /rr.2]  in 


ft  -  q2/2m  -  q  kFj  /m 


(3.  2.  27) 


These  principal  parts  are  evaluated  here  so  that  when  the  imaginary  parts  of 
these  integrals  vanish  (which  they  do  for  a  certain  range  of  the  q  vector 
specified  by  the  logarithmic  singularities  in  (3.  2.  27)),  they  directly  give  the 
information  concerning  the  real  poles  in  the  u-plane.  These  real  poles  in 
the  correlation  function  correspond  to  resonant  bound  states.  When  the 
imaginary  parts  of  these  integrals  are  finite,  the  reality  of  u  is  lost;  the 
bound  states  have  finite  lifetimes  in  these  approximations.  The  physical 
mechanism  for  these  lifetimes  has  already  been  outlined  in  Chapter  I,  and 
will  not  be  repeated  here. 
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For  small  q  ,  these  expressions  can  be  further  simplified  by  making 
approximate  evaluations: 

Ajj  /mu^  ;  Bjj  ~  -q^Nj  /mu^  .  (3.  2.  28) 

The  evaluation  of  Cjj  for  this  limit  is  a  little  more  subtle  because  of  the 

presence  of  2K0'?.  Writing  fl  =  (<j-2K0'£)  it  is  here  given  up  to  order 
2  —3 

(q  /SI  )  the  reason  for  which  becomes  evident  when  the  spin-wave  pole  is 


evaluated. 


(3.  2.  29) 


The  poles  of  y  ^  ,  7-j^  ,  and  are  *orm 

1  -  2/  2  Ajj  Bjj  =  0=  1  -  l/2q4N|N|  /m^u* 


(3.  2.  30) 


There  is  thus  a  sound  mode  in  the  medium  traveling  with  a  velocity  given  by 

F 

(vj_,  is  the  Fermi  velocity  k^,/m).  For  the  P  state,  KO'/e^.  -  2/3  ,  this 
is  Vp/  y/’S  which  is  the  usual  result.  It  must  be  remarked  here  that  in  the 
F  state,  Bjj  =  0  (only  t  states  are  occupied)  so  that  C  jj  has  no  zeros  or 
y  has  no  poles  showing  that  there  is  no  sound  mode  propagation.  However, 
the  spin  waves  are  seen  to  exist  in  this  state.  Note  that  a  plasma  mode  did 
exist  in  the  F  state  in  the  Coulomb  case,  in  (A)  .  This  result  will  be 


1 


range  interaction,  the  particles  interact  only. when  they  are  very  close  to  each 

other.  In  the  F  state  all  the  particles  have  up  spin  and  so  the  Pauli  principle 

comes  into  play  and  precludes  the  interactions  entirely.  (The  sound  mode  is 

4  found  in  yQ  ,  which,as  was  pointed  out  earlier,  may  be  thought  of  as  singlet 

state  of  the  pair;  in  the  present  instance,  no  singlet  state  can  be  formed. 

This  is  another  reason  for  not  finding  a  sound  mode  in  the  medium.)  The 

longitudinal  susceptibility  also  does  not  show  any  resonant  response  (it  did 

in  the  Coulomb  case  in  (A))  because  the  spins  are  already  lined  up  and  no 

longitudinal  field  can  fluctuate  it  any  further  because  the  interactions  are  of 

zero  range.  However,  the  transverse  susceptibility  shows  the  spin-wave 

mode  as  must  be  the  case  and  as  is  reasonable  otherwise.  Note  in  this  case, 

Y  ^  =  *  3I3  =  -  Ajj  which  can  be  expected  because  N  =  for, 

N  =  N|  +  Nj  =  N|  (F  state)  and  *3  =  N|  -  Nj  -  Nj  (F  state). 

The  spin-wave  mode  will  now  be  discussed  in  some  detail.  This 

requires  the  evaluation  of  the  zeros  of  (1  -  in  the  limit  of  small  q 

2  3 

and  small  u  ,  such  that  terms  of  order  (uq  )  and  (u  )  are  neglected.  Then 
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The  stability  of  the  spin-wave  mode  may  be  ascertained  by  imposing  the 


condition  that  its  frequency  be  positive.  This  criterion  gives 

K$' /eF  >  [(1  +  r)5/3  -  (1  -  if)5/3  ]/5?  . 

For  the  fully  ferromagnetic  case  this  is 

K0‘/e~  >  25/3/5 ~  0.  635 


(3.  2.  32) 


(3.  2.  33) 


From  Chapter  II,  the  ground  state  is  an  F  state  if 


K0'/ey  >  2'1/3  —  0. 


(3.2.34) 


and  the  ground  state  is  UF  if  K0'/cj.  >  2/3  0.  667.  Thus,  whenever  the 

ground  state  is  cither  UF  or  F  the  spin  waves  arc  stable  also.  This  was  first 
noticed  by  Herring  [23]  and  again  recently  by  Fukuda  [12]  . 

There  is  a  second  type  of  instability  of  spin  waves.  Here  the  spin 
waves  merge  with  the  scattering  states  thus  acquiring  a  finite  lifetime  even 
in  RPA.  This  happens  only  beyond  a  certain  maximum  wave  vector,  qmax 
of  the  spin  wave,  becaure  X  +  will  then  have  complex  poles.  The  imaginary 
part  of  begins  to  appear  first  for  qmax  and  given  by 


Ucu ,  -  q^  /2m  +  2K0'£  -  q  k  »  /r 
SW  1  max'  fci 


(3.  2.  35) 


Calling  q  =  Q  k  —  .  one  arrives  at  a  transcendental  equation  for 

6  nmax  max  F  ’ 

Q  ,  obtained  by  putting  (3.  2.  35)  in  the  spin  wave  dispersion  law 
1  =  77cu  (3.  2.  27): 
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=  |  (K0'AF) 


^(1  +  f)2/3+2(l-  f)1' 


-(i  +  o1/3)  + 


Q  -  (1 

max  ' 


+  o1/3-(i  -r)I/3  I] 
+  r)1/3+(i  -r)I//3  / 


(3.2.26) 


This  has  to  be  solved  numerically.  This  numerical  solution  is  displayed  in 
Fig.  6.  F rom  (3.  2.  26)  and  from  the  schematic  diagram  in  Fig.  6,  it  is  seen 
that  Qmax  has  an  upper  bound  at  ( 1  +  f) ^  ^  ( 1  -  ?)  *  ^ .  In  this  figure, 

F0'/eF  is  taken  to  be  an  independent  parameter  not  determined  by  the  ground 
state  conditions.  The  only  relevant  parameter  of  course  is  the  one  pertaining 
to  the  ground  state  and  the  corresponding  curve  is  shown  in  Fig.  6  as  a  broken- 
line  curve.  For  £  =  1,  the  F  state,  this  gives  an  exact  answer,  namely 


Qmax=  V2  (K0'Af)2- 


For  such  a  lx>  the  spin -wave  frequency  is 


(3.  2.  27) 


Ucw  (Q  k  )  -  f, r 

bw  max  £  F 


/K0_'\  _9_  2 1  /3  Kg' 

l  r)  ® 


(3.  2.  28) 


)»-  I 

The  spin-wave  stabil.ty  criterion  applied  to  this  shows  that  (Q  k_) 

SW  max  F 

is  stable  only  if  K g'/cp  >  0.71.  This  implies  that  for  stable  F  states  for 

which  KO'/fr  >  0.794  the  spin  waves  a  re  stable  up  to  Q  but  beyond 

*  max  7 

QJTiax  they  are  unstable,  being  scattered  into  the  individual  particle  excitations, 

even  though  the  infinitesimal  spin  wave  stability  criterion  permits  spin-waves 

to  pei  s  ist.  If  we  take  the  stable  F  state  for  whi'h  K0'/e„  =  2" 1 

7  F 

(3.  2.  28)  gives  the  spin  waves  up  to  Q  =  0.  945  and  its  energy  is 

max 


)  —  2 

o  W  a  x 


-/l 6  or  <  i  /16  . 


(K0*/«f) 


DOTTED-LINE  CURVE  REFERS  TO  K0'/*F  CONSISTENT  WITH 
THAT  EVALUATED  FOR  THE  GROUND  STATE  FOR  VARIOUS 


I 
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It  must  be  mentioned  that  using  an  approximate  method  for  evaluating 
2  \ 

the.  coefficient  of  q  for  (Jsw-  Herring  [23]  ,  Thompson  [24]  ,  and  Fukuda  [12] 
had  each  obtained  the  same  coefficient  as  given  here  for  the  F  state. 

Thompson  [24]  had  also  derived  a  value  of  Qmaae  for  the  F  state,  in  agreement 
with  that  given  here.  Baym  [27]  had  only  estimated  this  Qmax  for  f  <  <  1 
and  arrived  at  a  value  of  Qmax  =  2-1/3  (Kfl'/fj')?  >  which,  as  expected,  is 
smaller  than  that  for  the  F  state.  The  schematic  diagram  of  how  Q 

max 

comes  about  is  given  in  Fig.  la  as  well  as  Fig.  6  which  has  been  explained 
already.  In  this  connection,  an  aspect  of  finite  temperatures  must  be 
mentioned.  For  this  case,  the  integrals  A,  B,  C  all  acquire  finite  imaginary 
parts  and,  hence,  the  collective  modes  or,  in  particular,  the  spin-wave  mode, 
will  be  slightly  damped.  Physically  this  is  reasonable  for,  at  finite  temperatures, 
particle-particle  scatterings  occurring  into  empty  states  below  the  top  of  the 
Fermi  distribution  are  more  frequent  and  these  tend  to  destroy  the  coherence 
of  the  bound  state. 

In  the  light  of  these  results,  one  may  remark  on  the  nature  of  the 
difficulty  in  computing  the  total  magnetization  of  the  system.  We  have  two 
relations  for  calculating  the  magnetization: 

<  <r3  >  =  -  i  [tr  t-j  G  ( l  1  +)  ] 

and 

<<r3(l)>  6(3)(l-2)  =  <^[<r+(l),  <r_  (2)  ]_  > 
at  equal  times.  The  direct  evaluation  of  <  >  from  the  single  particle 

Green's  function  in  the  interacting  system  is  very  hard  and  so  may  be  used 
only  for  very  approximate  calculations  at  high  temperatures  where  the 


ARPA-1 1 


III-26 


interactions  do  not  contribute  to  spin  waves.  But  at  very  low  temperatures, 
the  correlations  are  important  and  one  may  use  the  second  relation  which  can 
be  computed  to  a  much  higher  accuracy  than  G  .  It  is  not  clear  at  the  present 
time  how  these  two  have  to  be  reconciled  to  arrive  at  the  correct  evaluation 
of  the  magnetization  valid  at  all  temperatures.  It  is  clearly  beset  with  a 
considerable  difficulty  of  self-consistency. 

In  the  next  section  we  will  discuss  the  existence  of  spin  waves  in  a 
Yukawa  gas.  This  does  not  follow  as  directly  from  the  equations  as  for  the 
zero-range  gas.  One  needs  to  resort  to  certain  transformations  of  the  given 
equations. 

C.  Yukawa  Gas 

Let  us  examine  the  Eq.  (3.  2.  1) 

r.  (  k  ;  q )  =  -  Ti  -  q  )  y.  (  q )  - 

?!-  ?!  GMq)ri(q;  q)  G(q)  a  (k;  q-k)  (3.2.  i«) 

\CTI)  b  (q  -  k) 

where  we  have  used  T  -  A  /Z  .  In  this  we  take  A  =  -1  so  that 

o  o  o  ^ 

£  =  &  RPA  ’  w'ien  evaluated  for  u  =  0  and  for  small  wave  vectors 

gives  ~  1  +  k^,p-/k  ;  k^.^.  is  the  usual  Thomas -Fermi  screening. 

Hence  W /Z  has  the  appearance  of  a  Yukawa  potential.  Thus,  the  Yukawa 
potential  can  be  used  consistently  with  the  ground-3tate  properties  in 
examining  the  collective  excitations.  In  case  (A),  if  exchange  is  taken 
into  account  by  an  extended  RPA  of  the  type  suggested  earlier,  it  shifts  the 
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pole  and  also  gives  a  finite  lifetime  to  the  collective  states.  This  is  as 
expected  since  taking  into  account  more  interactions  would  tend  to  degenerate 
the  coherence  of  the  bound  state.  The  use  of  "plasma"  RPA  with  a  Yukawa 
potential,  however,  gives  results  very  similar  to  those  in  case  (A)  -  a  plasma 
pole  or  a  zero  sound  mode  depending  on  whether  the  screening  is  small  or 
large, the  transition  being  gradual.  But  spin  waves  will  not  appear  in  this 
approximation.  A  more  careful  calculation  involving  the  exchange  contributions 
is  needed  to  derive  spin  waves.  The  spin-wave  problem  is  now  examined  here 
in  detail. 

To  study  the  spin  waves,  we  investigate  the  equation  for  and  only 
the  component  ( ) j  |  in  this  equation  contains  an  inhomogeneous  term,  since 
the  inhomogeneity  is  caused  by  the  external  field  and  we  are  seeking  resonant 
response  to  an  external  field.  Then  equation  (3.2.  L')  takes  the  form 

( 2/y  =  ^Yukawa1 

r2|  }  q)  =  '  1  + 

+  i  J 2/y  (k-q)  G  j  j  W+q)  r2|  |  (<i;  q)  G|  |  (q)  •  (3.2.39) 

To  study  the  collective  mode  we  have  to  consider  the  homogeneous  counterpart 
of  this  equation  which  we  will  now  examine.  This  is  equivalent  to  studying  the 
pole  structure  of  .  From  (3.  1.12)  if  the  usual  RPA  is  used,  with  r  (12;  3) 
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“  *h'  tWrd  ‘I”1  *°  -«(4,(l-2)-  S<4,(1-3)  and  neglecting  STJtt  U,  it  can  b. 

^..Uy  verified  that  1^(12;  3)  is  of  the  form  « (t,  -  t.,)  !>,•  (12;  3 )  .0  that  in 
Fourier  .pace,  one  may  ae.ert  that  in  (3.  2.  39),  r,<k;  q)  t.  independent 
of  k0  ,  the  frequency  part  of  k  .  Then  the  ^  integration  in  (3.  2.  39)  can 
be  carried  out,  which  is  just  C  as  in  (3.  2.  17c): 

r>!  1 1 Ik !  q  ’  *  I  &  *r  * •  5>  , 

[jJ  -  (J_(£+3)  +  u+(3>  _|  2»* 

(3.  2.  40) 

(homogeneous  part  of  (3.  2.  391)  where  us±(q)  are  evaluated  with  the  Yukawa 
potential  as  in  Chapter  II  (2.  5.  10).  It  is  interesting  to  note  that  in  the  limit 
"  ‘  O’  9  =  0  (in  this  order  of  limits),  (3.  2.  40)  becomes  the  ground  state 
condition,  namely  (2.2.  15).  Hence,  in  the  P  state  limit,  the  two  equation, 
should  yield  the  same  results.  We  now  make  an  algebraic  transformation 


Then 


r2JfO<;  q)  =  [u-cj  (k+3)  +  u+(k)J  r2jj(k;q) 

(u  -  (J_  (k  f  a)  +  u+  (k))  r2|j  (k  ;  q )  =  - 

'I  ("i,(Va»  -  nF  (u.  (3*3)))  r2||(q;  q) 


(3.2.41) 


(3.  2.  42) 


u+(k)  =  k  /2m  -  J  g/  (k  -  a)  n  (<[) 


(3.  2.  43) 


/ 
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F ukuda  [12]  arrived  at  an  equation  of  the  form  (3.  2.  42)  using  an  equation  of 

motion  method  for  the  spin  flip  operator  in  the  RPA  which  involves  the 

exchange  term  and  he  worked  out  the  Coulomb  case  in  the  F  state.  The 

present  approach  follows  a  procedure  somewhat  similar  to  that  of  Fukuda  in 

solving  (3.  2.  42)  for  a  Yukawa  potential,  but  for  all  f  and  f  . 

We  will  first  indicate  how  (3.  2.  42)  can  be  solved  for  quite  general 

and  then  specialize  to  the  Yukawa  potential.  Let  y  =  k^c F ,  y  =  q/k  , 

2 


the 


v  =  'lAj.  ,  -  =  (kf  /2m)  v  and  or^  =  mc2/^^.  and  also  cliooae  x  as 
z  axis  of  a  reference  coordinate  system.  Then,  since  cos  9  =  )*  /2  YlQ(0^  ), 

being  the  usual  spherical  harmonics,  we  have 


v  -  x 


-2xy(-j-)  f  Y1()(y)  +  (  -- ;£)  y  d3y  V  ( |y-y  J  )[nF(u_(x+y))  nF(u+(y>ll^j|(y;x) 


ar  f  3_ 

=  ( )  J  d  y  v  ( I y-y  I )  [nF(u_  (*+y))-  nF  (u+  (7))j  r2j  |  (7;  x) 


(3.  2.  44) 


where  we  have  written  for  convenience,  *1?  (y  y)  =  - n  °L  V  ( |  y-y  | >  .  Let 

kF 

x  -  0  in  thi3.  Multiplying  the  resulting  expression  through  by 

(nF(w_(y))  -  nF  (u+(y))  )  and  integrating  over  y  ,  we  see  that  a  solution  exists 
with 

v  =  0  and  f  j  j  (y,  0)  =  a  constant  .  (3.  2.  45) 

Multiplying  (3.  2.  44)  through  by  (nF(u_(y  +  x))  -  n^  (tJ+(y)))and  integrating  over 
y  we  get 


ARPA-ll 


UI-30 


y  n  .3  _  \ 

(v-xZ)  j  — - £j-  (nF(u_(x+y)-nF(u+(y))  )r2jj(yjx) 

4| 

-2x  (-^-)  •'‘SlSr  yY10(9>  (nF(w>(x+y))-nF(w+(y)))r2J|  (y;x)  + 

(2l) 

+  ( — 2~)  d  y-  V(|y-y|)  (nF(u_(x+y))-nF(u^(y))}(nF(u_(x+y|-nF(fc^(y))) 


r?j  j(y;  x) 


=  { — 2“)^*^  ^  y  ^  y  V(  |y-y  |  )(nF(<j_(x+y))-nF(iJh(y)))(nF(u_(x+y®-nF(u+(y))) 

r?.l  t (y;  x> 

Interchanging  y  and  y  in  the  RHS,  since  V  is  in  general  symmetric,  this 
term  cancels  the  last  term  on  the  LHS  and  hence  one  arrives  at  a  general 
expression  for  v  for  all  x  ,  lj  ,  and  $  : 


^  3 

+  2x  (-^J-)1/2  ^  y  Y  io(y)fnF(u-(x+y)>*nF{u+(y))1  r2|  t  {yi  X) 


\  - \  (n (j_(x  +  y))-nF(u+(y))J  I^j  j  (y;  x) 

(3.2.46) 

To  find  the  coefficient  of  x2,  we  expand  (nF(u_(x+y))-nF(u+(y))  )  and 
I\,j  |  (yj  x)  in  powers  of  x  taking  note  of  (3.  2.  45) 

OO  OO  ''v 

nF((j  (x+y))-nF(u^(y))^n(0)(y)+^  ^  x‘  Y^Q(y)  n^^(y) 

i=  1  f^O 


r2l  |  *>  *  rzl'|0)+  II  J  r2l  t 

j=l  fm' 


(3.2.47) 


Ill- 31 


ARPA-ll 


Yjg,,  appears  in  the  first  expression  as  the  left-iiand  side  contains 
explicitly  only  x»  y  (or  ^ )  while  this  is  not  so  clear  in  the  second  expression. 
Working  then  only  to  lowest  order,  we  find 


,,/z  I y3"(0)(y)r2['|)i0(y)dy ^ I dy 

L  4¥r2i|0)  J  y2n^0)(y)dy  / 


(3.  2.  48) 


Now  clearly  at  T  =  0°K, 

nF(ui_(x+y»  -  nF(u+(y))  =  [nF(u Jy))  -  np  (u+ (y))J 


-  2  (  ^  )1/2  y  Yln$)  6((l-C)2/3-y2)  f  .. 


so  that 


n  vy)  =  nF(u.  (y))  -  nF(w+  (y)) 


nj(1)(y)  -  -  2  (-^-)1/2  y  6p-K)Z/3  -  y2  ) 


(3.  2.49) 


Choosing 


r2||(0>  "  <3/4tt  )1/<S  (ars/ir) 


(3.  2.  50) 


we  finally  get 


usw  =  <1  /z™?  V  -  (~)  S  <?,?> 


with 


(3.  2.  51) 


S  (?,  S)  =  y  dy  y3  (nF(U_(y))-nF(u+(y))  )  ^(y) 
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Now  our  task  is  to  find  j  ( y)  .  This  is  done  by  using  {3.  2.  47)  in  (3.  2.44) 
after  multiplying  it  by  [n^lu  (x+y$-nj.{u+(y))]  ,  along  with  the  expansion 


V  =  X  +  .  .  . 


(3.  2.  52) 


and  then  equating  the  coefficients  of  like  powers  of  x  on  both  sides  of  the 
resulting  expression.  Before  doing  this,  let  us  note  that  formally 


v  (|y-yl)  =  —7  y  Q|..  (y;y)  ( 


(3.  2.  53) 


+1 

—  C  2 

where  p  ‘is  the  angle  between  y  and  y  .  Using  j  P^{p)  P|,(p)  dp=£f£-j  V 
we  find 


r  i 

4-Q|< y;y)=y  p/p)  v([y-y|) 


(3.  2.  53a) 


In  (3.  2.  53),  referring  y  and  y  to  a  system  where  x  is  the  z  axis, 


P/«  (p)  4lT  ^  *  2i  +  1  ^  Yf"m"  Vm"(y) 


so  that  (3.  2.  53)  takes  the  form 


v  (  I y-y i )  =  — z~y  Qi..  (y;y)  Y|i.m.i(y)  • 

yy  ,Vm» 

Here  Y,  'a  are  such  that 

tm 


(3.  2.  54) 
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^  <*>  v*  m,  (?> = ( 
+1  2* 

where  J  dy  s  §  d#i  as  is  usually  taken.  In  the  Yukawa  case, 

(3.  2.  53a)  takes  a  particularly  simple  form 


(3.  2.  55) 


1  _  r  p*(m) 

— r  Q#(y;y)  =  \ <4i  — ? — y 

YY  1  fvZ+7Z+ 


But  since  ( t )  = 


_i  [y  +y  +  5  - 2yy /*] 


function  of  the  second  kind 


if,  I 
5  T  J  "*  TT 


P» 


wc  find 


Q^yjy)  =  Qf 


2,_2X  If2 

r  +  y  +  f 


2  y  y 


(3.  2.  56) 


The  expansion  (3.  2.  53)  was  anticipated  in  view  of  this.  Using  these,  we 
get  for  the  coefficient  of  x 

-2y  ( 3-1 1  /2  Y10<?)  r2|,(0)+  2  (-jy )  £  v(,m,<r>  r2|(,  (y)f  y  ->(0)(y)Qo(y;  y)dy 


=  3  <?F>  X  YiW<Y>  i  r  ",0>(y>  o,.(yiv>ra|V/>-m. 


(3.  2.  57) 


*  ,A 


Multiplying  this  through  by  (y)  and  integrating  over  y  ,  we  get 


i 


ARPA-1 1 


III- 34 


-2y(4>l/2r2|  4mot21^>lj  yn'“'(7)Q0(y;7)dyj  rrfj£iW  • 

=  2<^>  fy  »<°>  (y) Q|(y;y) (y) iy  .  t3-2-58' 

_  (1) 

Note  that  for  !  /  1  ,  m  /  O  ,  this  gives  homogeneous  equations  for  ^  which 

can,  therefore,  be  taken  to  be  zero  without  loss  of  generality.  Also  only 
r2|  |(^  (y)  appears  in  the  expression  (3.2.  51)  for  Ugyp  Thus,  we  consider  the 
equation  for  f  =  1  m  =  0  only.  Using  (3.  2.  50)  for  j0)  ,  we  finally  have 


l^?>  ft  7"(0)(7)Qo<y;7Myj  r2,f<”  W 


-y2  + 


|  j  7  n<0)  (7)  Q0(y,-y)df  |  r2|||(0'>(y>  - 

=  ^  7  n(0)  (7)  Qi  (y;y)  ^2\  fVo  ^  d*  * 


(3.  2.  59) 


This  equation  must  be  solved  for  I^j  |  to  arrive  at  an  explicit  exprcssion- 
for  the  coefficient  of  q2  in  Note  that  the  above  framework  is  so  far 

independent  of  the  choice  of  interaction,  and  hence,  is  a  proof  of  the  existence 
of  spin  waves  for  quite  general  interactions. 

Now  we  will  specialize  V  to  be  a  Yukawa  potential,  so  that  (3.  2.  56) 
has  to  be  used  in  (3.  2.  59).  Noting  in  this  case  that 


Qt(y;y)  = 


2  -2  ,  e2 
'  +  V  +  ? 


4  y  y 


•2  +  (y  +  y)2 


r  +  (y-y>‘ 


(3.  2.  60) 
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and  since  the  integral  on  the  left-hand  side  can  be  explicitly  done,  we  finally 


have 


3  2 

y  +Y 


■) 


r+(y-d +c>  /Jr 


g2t(»-r)2/3-y2  L/^«(y(»-»‘^V  \ 
* I  \  c'My-O-iwV  / 


» t  [u.-1  | 


+  tan 


+  tan 


-1 


onL/,-y 


d-n'^-y 


1 

I) 


ilflo  <*> 


•s 


dy 


2  —2 

L±Zl±X— 

4 


tn  I  IJ&tK 


-  y  y 


(nF(o+(y)}-nF(w  (y>)) 


r  ..(D 


2j  1 10  (y) 


(3.  2.  61) 


Now  note  that  if  /l0  =  -  ^j^^ty),  this  equation  still  holds:  The 

equation  (3.  2.  61)  is  an  inhomogeneous  Fredholm  equation  of  the  second  kind 
and  hence  admits  of  a  power  series  solution.  Furthermore,  this  has  an 
iterative  solution,  the  first  iterate  being  the  inhomogeneous  part.  Since 
the  coefficient  of  ( y )  the  left-hand  side  of  (3.  2.  61)  can  be  shown 

to  be  positive  for  positive  y  ,  it  is  clear  that  the  solution  q  (y)  has  all 


l 
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its  terms  negative.  In  view  of  the  symmetry  noted  earlier,  we  may  seek  a 
power  series  solution  in  the  form 


-  (1) 
r2ltl0 


(y)  -  f 


'2n+l 


(3.2.62) 


The  integrations  on  the  right-hand  aide  can  all  be  done  in  principle  and  thus  a 
linear  expression  Involving  the  coefficients  ^  an^  obtained.  Equating  the 
like  powers  of  y  on  either  side  of  this  identity,  we  get  a  set  of  simultaneous 
linear  inhomogeneous  algebraic  expressions  for  {an^  »  which  can,  therefore, 
be  determined  in  principle.  Once  these  coefficients  are  determined, 
is  obtained  and  hence,  S  ($,  ?)  .  In  view  of  the  cumbersome  nature  of  this 
equation  (3.  2.  61),  we  have,  in  the  case  of  arbitrary  magnetization,  chosen 
a  single  term  expansion  in  (3.  2.  62).  Even  this  gives  a  very  complicated 
expression  for  S  (£,  £)■  We  here  give  S  in  the  Coulomb  and  Stoner  limits  for 

the  sake  of  simplicity  ; 


USR  =  a_J  i  _ 

u  SW  2m?  ) 


(IK)5/3  -  0-g)5/3 

5? 


(3.  2.  63a) 


whe  re 


_  (1) 
r2lt  10 


Y  for  C  -  oo 


c  s2  r,  fun)5/3  -  (i-g)^- 

USW  “  *Zm?  S  10  ar3  {l  +  ?)l/3  .  (l-?)  ^ 


where 


r2HlO(y)--  2  ((1  +  ?)1/3  -  (l-O1/3! 


for  ?  =  0 


(3.  2.  63b) 
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Note  that  (3.  2.  63a)  is  precisely  the  same  as  that  obtained  earlier  by  a  rigorous 
method  not  involving  the  steps  suggested  above.  In  the  paramagnetic  limit,  using 

i 

the  equilibrium  K0’/*p  and  ar,»  respectively,  it  is  found  that  the  spin  waves 

are  unstable  as  must  indeed  be  the  case.  Again  for  the  F  case,  (?-l), 

in  both  cases  are  positive  when  the  corresponding  equilibrium  coupling  strengths 

are  used  as  they  should  be. 

For  the  ferromagnetic  case,  however,  we  have  carried  out  a  three-term 
expansion  in  (3.  2.  62)  and  again  a  very  cumbersome  expression  for  S(?)  for 
all  £'s  obtained.  For  the  sake  of  simplicity  we  again  give  the  Gimkuub  and 
the  Bloch  limits  only.  In  the  Stoner  case,  a 3  and  a&  are  identically  zero 
to  [0(£2)]  ,  as  this  is  the  only  meaningful  term  here.  In  the  Coulomb  case, 


however,  evaluating  (a^a^  gives  S  ($--0}  ?=1)  =  0.3517  and  hence 

c(r=n  2  /,  5.344  , 

•Jew  =  <1  /2m  ' 1 - r - >  ‘ 


(3.  2.  64) 


This  must  be  compared  with  Herring's  result 

c(r=l)  2  /,  ,  5.485  . 

USW  =  q  /2m  (1  p— - 


Fukuda  [12]  gives  5  136  for  this  and  this  is  wrong  since  he  made  an  error 
in  computing  a&  .  It  may  be  remarked  that  a  two-term  expansion  gives 
5.  166  in  agreement  with  Fukuda.  It  should  now  be  stated  that  using  Herring  s 
spin-wave  stability  criterion  using  the  equilibrium  value  rg  =  5.  45  (Bloch[3]) 
we  find  that  our  calculation  gives  a  positive  u)gW  whereas  Herring's  gives  an 
unstable  spin  wave  (negative  UgW).  In  view  of  this,  we  feel  that  our 


calculation  is  better  than  Herring's. 
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In  the  next  section,  the  case  of  Bloch  electrons  will  be  dealt  with  in 


some  detail. 


3.  The  Case  of  Bloch  Electrons 


In  Section  4  of  Chapter  II,  suitable  modifications  for  treating  the  single 
particle  states  in  a  periodic  potential  were  outlined.  In  conjunction  with  . 

(2.  4.  2),  let  us  define  the  vertex  parts  in  the  Bloch  scheme  as  follows: 

r.{/fk;q)=y  b^+q(j)  bilk(l')  eiq'2  r.(ll';  2)  d4ld4l*  d42  (3.3.1) 

IBZ 

f:  band  index.  Then  the  equations  in  the  RPA  including  exchange  can  be 


written 


r.{/f'k;q)  =  -  Ti  <  fkf qjf'k  >  - 

n/k+q  f  q  -  4— 

^f.q+qfk  [tr{G  (q+q)r(V2q;q)G  (q)}]  ^ 

4  J  *  *  £  \  *-7T ) 


+  iT  f  ^^?i2+q-G^(q+^ri(il/2q;q) 
ft  J  (2ir)4  ikfiq+q  1  2 


(3.  3.  2) 


Here 


v‘‘k  Si =  b*1k^i'b’V2<'-2»b,3k3<ri>v4(ri,,',Tri)',3r> d3r' 


(3.  3.  3) 
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i.k.  i.k, 

Using  the  periodic  properties  of  b-.  ,  it  is  seen  that  V  /.  1  is  non-zero 

only  if  k|  +_k£  =  ^3  ^  JI4  (mo<l  K)  (K:  reciprocal  lattice  vector).  Neglecting 
Umklapp  processes  (see  Adler  [46]  for  their  effects)  and  using  a  Fourier 
representation  ioxDir^-r^)  one  arrives  at 


^ iV  =  (2r)3<l1k1|  f4k4  >  <  iz k2 1  f3k3  >  2/tkj -  k4) 

a(3)(k1  +k2  -  k3  -k4) 

Defining 

T{  (ff'kjq)  =  </k+q|i'k  >  ^'(ff'kjq) 
and  U3ing  (3.  3.  4)  these  equations  can  be  recast  into  the  form 


(3  3  4) 


(3.3.  5) 


r.'ff/'kjq)  =  -  ri  -  Vfq)^3^)  + 


+  >T  f  Vlk-ilfc/lU'kq;  »  (,q)  G(  (q+5)r  '(<,(2q;  q)  G. 

L-J  u  (27r)  *1  *  c  lz 


*1*2 


(3.  3.  6) 


where 


<  £k+q[f  q  +  q>  «  q+q[f  qXf  q  [  <'k  > 

kqjfjf^r  <lk  +  q|fk  > - (3.3.7) 


and 


YiB(q)  "  f  d  ^  4  \<  *tq+q  l*?q  >1 
1  u  J  Uv)  2 


V2 


[tr  [C^tq+q)1’.’  «ll2q;q)  G^q)]  ] 


(3.  3.6) 


Finally 
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*33  ^ 


+ 1  ^(q).Tr3B{q) 


£  J  |<  iiq+5|l2q  >|2  [fcr^G^tq+qjA^lj^q;  q)  G^lq)]  ] 


(3.  3.  13) 


X  ®<q)  =  i£  J  -^V  |<liq+q|l2q>|2  [tr[rf  G^Cq+qjA^lj^qiq)  ]  1  + 


+  i^(q)T1B(q)^  V  ^l< 
l\  lz 


^q+ql  £2q>i  2  [tr  ^  (q+q)A0(<i<2q;  q)G|^(q)  j  ] 


(3.  3.  14) 


In  the  'plasma'  RPA 


B  ,  „  B. 


CB(q)  =  1  +  V(q)  (Aj  +  Bj  ) 


=  1  +  i  V(q)  y  C  ^\-|<*,q+5l<2C>|2('r{G  (qt5)G  (q)}  ] 

J  I7v\  1  *• 


(3.  3.  15) 
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and  similarly  other  expressions  for  X33  ,  X+_  ,  etc.  One  of  the  important 
features  of  the  Bloch  case  is  the  appearance  of  additional  sums  on  the  band 
indices  and  the  overlap  integrals.  The  result  for  CB(q)  in  the  unpolarizcd 
case  is  that  first  derived  by  Ehrenreich  and  Cohen  [40]  using  the  equation  of 
motion  method.  These  results  here  are  generalized  to  include  the  polarizations. 
The  above  results  are  valid  only  when  the  local  field  corrections  can  be  neglected 
since  the  Umklapp  processes  have  been  neglected  entirely  (Adler  [44]).  This 


neglect  of  local  field 


CoiTtclionH  is  justified  when  the  bands  are  broad  and  far 


apart.  In  Appendix  B  the  relation  between  these  equations  and  those  that  can 
be  derived  by  the  equation  of  motion  method  will  be  given. 


In  view  of  the  complexity  of  the  results  derived  above,  only  the 
existence  of  spin  waves  will  be  discussed  for  the  case  of  a  single  band,  not 
overlapping  with  any  others.  As  in  Section  2  of  this  chapter,  the  equation 
describing  the  spin  flip  oscillation  in  the  RPA  including  exchange  is  written 
in  the  form  (3.  2.  42)  for  the  F  state  (we  have  treated  only  F  state  here  to 


keep  the  complications  to  a  minimum): 


(u  -  ^k(q))  J  (k;  q)  =  -  J  ^  (k-q'J'Zi/fkq;  q)  j  (q;  q) 

ek(q)  =  dh+q)  -  e(k)  -  f  --ij-  (<k|q  >|2  1/(k-^)  (3.3.16) 

u  <2tr) 


where  e(k)  is  the  band  energy. 
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In  the  almost  free  electron  limit,  <  k|q  >  ~  1,  and  so^tS  ~  If  since 
2  2 

4/  (k)  =  4ne  /k  the  arguments  of-ihe^ previous  section,  case  (C),  go  through 
and  hence,  spin  waves  exist  in  this  case. 

In  the  tight  binding  limit,  where  the  overlap  integrals  can  be  expressed 
in  terms  of  Wannier  orbitals,  a(£-_l  ),  is  the  lattice  site,  and  if  both  k  and 
are  near  zero, 


Ck|q>  =  ^ 


x  )  d^r 


unit  ceil 


■  Z  •  i_  1  J  *  te-i. 


unit  cell 


)  a  (r  -  f2)  d  r 


The  integral 


I  *V 


fj  }  a(r-f2)  d  r 


unit  cell 


is  the  overlap  between  two  Wannier  states,  Ijp  .  In  the  extreme  tigM  binding 


limit,  5  i  i ,  ,  and  if  k  and  k'  a  re 


••X-J 


d  1  so  that 


<  k|  q  >  =  O  ( |k  -q  |).  The  arguments  of  Section  2,  case  (C)  of  this  chapter 
are  applicable  so  that  the  spin  waves  are  seen  to  exist  in  this  case  also, 
since  (kq  q)~02(|k-q|)  and  the  equation  for  u  is 


(u-  ek(q))  r2J  }  (kj  q)  =  ^(k-q)  O2  (|k-q|)  T2j  j  (q;  q) 
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e.  (q)  =  e(k+q)  -  e(k)  -  f  A_i_  ^(k-q)  O2  (|k-q|)  . 

k  J  ( 2»r 

There  is  another  approximation  to  the  integral  (3.  3.  4)  suggested  by 
Antonoff  [29].  The  corresponding  results  re  here  derived  using  the  present 
technique,  since  they  follow  quite  easily  from  the  expressions  already  derived. 
Antonoff  makes  the  "degenerate  kernel’'  approximation  which  consists  in 
taking,  in  the  single  band  scheme 


k.  k? 
nt  1  2 

°  k  k 

k3  4 


Vex(k,-k3) 


(3.  3.  17) 


k.  k- 

What  is  done  here  is  to  write  .  in  terms  of  Wannier  orbitals,  and  in  the 

k3  4 

resulting  expression  retain  only  two-center  integrals.  Then  (3.  3.  2)  takes 
the  form  (in  the  one— band  case) 

r  (k;  q)  .  -  Ti  -  i  C  A-  (Vc(<,)  +  Vex(k-5  ) 

1  u  (27T) 

[tr  (q+q)  F  (q;  q)  G  (q)^  ]  + 

+  i  f  ?  i  (ty  (k-q)  +  *Z/  (q))  G  (q+q)  T  (q;  q)  G  (q)  (3.  3.  18) 

J  (2jr  )  c 

In  this  one  further  neglects  ^ex(k-^  in  the  second  term  and  2/.(k-q)  in 
the  third.  Then  a  calculation  similar  to  case  (B)  of  the  present  chapter 


follows : 
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=  -  2  Atti/(1  +  (2-*/ (q)  -  V  (q))  Am) 


(3.  3.  26) 


III  III  rt 

y2  =  T3  =o 


(3.  3.  27) 


X33m=-2Am/(1-,J,exWAIU>=2X“ 


(3.3.28) 


From  these  results  one  may  notice  the  following  features,  even  though  many 
terms  are  neglected  in  arriving  at  these  results.  Tlie  Coulomb  and  the  exchange 

i  —  —  1  —  .  n—  „  •>  —  ...  —  11  f  U  .  o  ^  w  ••  A 1  —  J M  m  f  1 1 v.  „  3  J  —  n  .  ..  —  An  3  3/  ^  «  ».V  «  — »  «  —  «1  ..  O  0 

VM  A  4t4  U4*  k»»v  V  V/  *  A  <-  44  Vii  4U31CV1V  tiw  #V  •  W  44  W  A  C  OilA  Jl 

•  *  C  X 

appears.  This  is  significant,  and  the  reason  why  in  the  "plasma"  RPA  X  + 
did  not  show  the  pole  structure  directly  and  why  it  did  in  case  (B)  now  appears 
obvious,  even  though  this  fact  that  exchange  has  to  be  taken  into  account  fully 
was  stressed  earlier.  In  case  (B), however,  we  had  taken  -  *2/  =  , 

which  when  substituted  in  the  above  set  of  expressions  give  back  those  found 
already.  The  results  of  case  (A)  follow  if  ^ex  *3  neglected  altogether. 

It  must  be  mentioned  here  that  if  the  "degenerate  kernel"  approximation 
is  made  in  the  Hamiltonian  itself  after  rewriting  it  in  the  Bloch  representation: 

H  =  Jj  *  (k)  akt  ak<r  + 


'till  <\ 

klk2k3k4ffl<r2 


(kl  +  k2=  k3  +k4) 


1  2  a  +  a  *  a 

3k4  Vl  V2  k3'2  Vl 


(3.3.29) 
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one  may  then  combine  the  terms  in  such  a  way  as  to  get  an  explicit  Heisenberg-, 
type  interaction.  Actually,  4^x(q)  goes  here  with  terms  which  can  be  cast 
into  the  form  <r(q).  <r (-q)  which  is  just  the  Heisenberg  result.  Now  if  RPA 
is  used  for  constructing  the  various  correlation  functions,  the  results  obtained 
here  a.re  recovered.  This  has  been  verified  by  the  author.  Antonoff  [29]  in  his 
thesis  neglected  ‘S/.tq)  entirely  and  obtained  spin  waves.  As  remarked  above, 
this  is  equivalent  to  assuming  a  Heisenberg  interaction  between  spins.  In  later 
work  with  F.  Engle rt  [29]  Antonoff  has  included  4/  (q)  and  obtained  the  results 
given  above.  They  used  an  equation  of  motion  method. 

In  summary,  the  possibility  of  spin-wave  excitations  in  an  itinerant 
model  is  here  shown  under  various  model  interactions  within  the  context  of 
RPA.  Moreover,  the  case  of  an  electron  gas  interacting  with  a  Yukawa 
potential  is  studied  explicitly  and  the  coefficient  of  q2  in  the  ferromagnetic 
state  is  determined.  The  results  are  in  good  agreement  with  those  obtained 
b'y  Herring  [22,  23]  in  the  appropriate  limits.  The  collective  excitations  for 
the  Overhauser  situation  are  not  examined  here. 
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There  are  many  objections  to  the  use  of  Yukawa  potential  (even  though 
potentials  of  this  type  with  Thomas-Fermi  screening  ar«f. suggestive)  in  such 

I 

calculations  as  are  undertaken  in  the  present  work.  The  first  major  objection 

is  that  it  is  an  arbitrary  choice  and  in  reality  the  situation  may  be  more  complex. 

;  * 

Secondly,  even  if  the  choice  be  arbitrary,  no  method  is  here  outlined  to 
determine  the  screening  parameter  from  first  principles.  Thirdly,  the 
screening  of  the  interaction  in  this  way  is  known  to  reduce  the  contribution  of 
the  exchange  energy  to  the  cohesive  energy  enormously  if  the  screening  radii 
are  chosen  to  give  agreement  with  the  observed  specific  heat  (£kp/'-'  10  cm  ). 
(See  Pines  [38].)  In  spite  of  these  objections  the  effect  of  screening  on  the 
various  properties  of  the  electron  gas  Is  most  easily  Investigated  analytically, 
by  using  a  Yukawa  potential. 

There  are  several  possible  directions  of  extension  of  the  present  work. 

Concrete  results  are  derived  only  for  T  =  0°K.  These  must  be  generalized 

to  finite  temperatures.  This  extension  should  settle  the  query  whethe  the 

3/2  2 

low  temperature  magnetization  is  the  Bloch  T  or  the  Stoner  T  It  is 

almost  clear  from  the  present  work,  that  at  very  low  temperatures  one  must 
have  T  '  behavior.  This  view  has  been  put  forward  by  Brooks  [6b]  and 
very  recently  by  Mattis  [47]  also.  There  have  been  3ome  diagrammatic 
analyses  of  this  problem  which  are  not  quite  rigorous  or  conclusive.  The 
equations  describing  the  magnetization  are  quite  nonlinear.  Some  self- 
consistent  way  of  calculating  the  magnetization  at  low  temperatures  must  be 


E-l 
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developed  as  was  stressed  at  the  end  of  the'last  chapter.  The  notion  often 
expressed  and  notably  by  Edwards  [25]  that  the  Stoner  and  Bloch  terms 
for  magnetization  and  specific  heat  will  both  be  present  is  almost  certainly 
,  wrong,  as  the  existence  of  qmax  for  spin  waves  clearly  indicates.  Even  at 
T  =  0  K  ,  the  problem  has  to  be  done  more  completely  by  including  all  the 
possible  processes  in  computing  the  ground-state  energy  of  the  system  as 


a  function  of  magnetization.  The  present  work  indicates  very  clearly  that  a 

mere  comparison  of  F  and  P  state  energies  is  certainly  not  sufficient  to 

determine  the  nature  of  the  ground  state,  as  was  done  by  Cooper  [13]  and 
Shimuzu  [11].  This  is  hard  only  in  the  sense  of  numerical  analysis  but  the 
way  U  can  be  handled  is  outlined  in  this  work,  for  the  expression  for  the 
dielectric  constant  is  here  derived  from  which  the  ground-state  energy  can  be 
computed  as  was  done  by  Cooper.  The  third  extension,  which  is  not 
considered  here  at  all,  is  the  collective  excitations  of  the  system  when 
Overhauser  type  ground  states  exist.  These  must  show  some  new  characteristics. 


Professor  P.  C.  Martin  considers  SOW 


as  an  itinerant  version  of  the  usual 


antiferromagnetism.  If  so,  one  may  expect  antiferromagnetic  spin-wave 
type  oscillations  in  the  SDW  case.  We  have  now  succeeded  in  solvine  this 

O 

problem  in  an  approximate  way;  and  the  results  on  this  and  other  matters 
concerning  itinerant  antiferromagnetism  will  be  dealt  with  in  a  separate 
technical  report  The  fourth  extension  is  a  more  thorough  investigation  of 
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the  results  given  here  for  Bloch  electrons.  The  fifth  problem  is  the  examination 
of  the  Overhauser  problem  once  again  in  the  Bloch  scheme.  The  investigation 
of  SDW  for  Bloch  waves  must  be  of  interest  in  real  solids  as  extrapolation  from 
the  free  gas  results  to  the  real  cases  is  often  misleading.  The  outstanding 
problem  of  all,  of  course,  is  a  complete  investigation,  even  within  HF,  of  all 
the  possible  ground  states  including  those  of  broken  symmetry.  The  phase 
diagrams  given  in  Figs.  2a,  b  do  not  include  the  superconducting  and  the  SDW 
ground  states  and  if  these  could  be  somehow  incorporated  in  this  plot  one  would 
have  a  more  complete  account  of  the  electron  gas,  even  at  T  =  0°K. 
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APPENDIX  A 

\  SPECTRAL  REPRESENTATION  OF  THE  CORRELATION  FUNCTIONS 

By  virtue  of  the  definition,  commutation  rules,  the  Hermiticity  condition, 
and  the  trace  structure  of  the  correlation  functions  defined  by  (3.  i.  2),  a  set  of 
formal  relations  is  obtained  for  their  spectral  weight  functions.  These  are 
given  here  without  derivation  (c.f.  Chapter  II,  section  1). 

App  ^rl  r2  *  u  ^  »  A33  ^  rl  t2  »  an<*  ^+-^rlr2*  u')  correspond  to  y  , 

Xjj f  snd  • 

App  (rir2;u)  =  App  (r2ri; 


I  z?  App,ri r2;">*o 


Lt 


App  (  rl_r2;  u)  =  1  c-0  *J_r2;  u+  ie)  '  yo  {  rl  r2;  u  _ie)l 

r  App{r 


<p  (2)  p(  1 )  > 


I,.,*  =f 

2  1  J 


"2i r  ,  pu 

(e r  - 


t0(12;  z)  =  j 


diJ  APP<ririiU’ 


(z  is  a  complex  variable) 


(A-  1) 
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A  33  r2  »  =  A33  <  r2  ri  »  *«• 

J  5TA33(rl  r2<  “>  -  0 


A33,rl  r2i  “I  *  *  “0  (X33(rl  r2i“+i  e  >  -  X  33<rl  ''2;  u-U  )) 
<J3(2)  t (1»|  .  f  *>.  A33<rjr2i“> 

3  *2  *1  J  21 


(J  -  t 


A+“  (rir2*u)  =  A-  +  (r2rli  ~u) 

I  7T  A+-  (ril  r2»'w>  =  <  *3<rji>  >  «(3)(rrr2) 

(equal  time  commutation  rule) 

[<r+(l),  <rJ2)]_  =  <r3(l) 

A+  (rIr2_'u)  =  le-0  [  X  +  .fr^r^u+ie)  -  X  +  .f^r^u  .  ie)] 


<  a^(Z)vJl) 


>|  ,  =  r_du  A+-(rl*2;^ 

Vli+  "  J  ^  p^.ij 


(A -2) 


U  -  z 


(A-3) 
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TUB  RELATION  BETWEEN  THE  EQUATION  OF  MOTION  METHOD  AND  THE 

GREEN'S  FUNCTION  METHOD 


Here  the  equation.  describing  the  variou.  collec.Iva  excitation.  given 
in  Chapter  IH,  will  be  rederived  by  extending  the  equation  of' motion  method  to 
polarized  medium  for  Bloch  electrons.  Incidentally, this  extend.  Wolff's[31] 
method  to  Bloch  electron,  and  al.o  include,  the  free  electron  re.ult  a.  a 
special  case.  The  Hamiltonian  i.  written  in  the  Bloch  representation: 


H=X  T  r#  a.  , 

f,k,  ^  1  >‘l  klVl 


t'j  y  y  y  +  + 

,lkll2k2  VlVi  'la2  ^ k4<4  k;<1<r‘  a''2'2''2  * 

(l<l+k2=k3+k4  (mod  K)) 

%2  |k;k2'ilW2>f<v  w  * 


U-'ql,+  'ak!'1'ak2M.'2|  tak2',iV<ll'2l  ’  U»,q‘) 

^Vl1  “Wc*  ""Vi1  akz  +  <ll/2l  >  LI3  (qj^T  ( 

\  .  A  J 

a  and  a  are  the  creation,  annihilation  operators  for  Bloch  electrons 
satisfying  the  usual  anticommutation  rules.  Then  the  equations  of  motion 
m  the  RPA  for  various  operators  are  constructed,  after  usin? 


(B-  1) 


$ 
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(2’,,3®(3,(kl+k2-k3-k4)«'H'l-k4) 


<k,i1lk4i4><k2l2lkJ<3>. 

{Umklapp  processes  are  neglected  as  in  the  text) 

r”‘  (k0'-  *•>  =  <*UV  *k+qto'  1  ~ 
f(J  -  u/oc+q)  +  ul.  (k)]  ^  (klf'  ■  q)  = 


=  (nF(u|+(k+q))  -  nF{(j^,  (k))]  <  k+ql|kl'  >  J 


-  uo(q)  -  u3(q) 


■'2/(q)  I  Z  1  <  >r"(q.i|f2;<1>  ♦ 


q'  v2  » 


Y  V  l/lk-q'MkH'qjq'i  f)  .. 

V  /  - : - r'r  (q'f.f  • 

£  v2  <q'*q'|h''2  >  12 


q>* 


is  as  defined  in  {3.  3.  7) 
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[u  -  uf"  (k  +  q)  +  U/("  (k)]  r*  *  <k W,  q)  = 


r 


3  lnF(W|‘(k  +  q))-nF(ufr  (k))]  <  kiq/|kf'  >  UQ(q)  +  (q) 

'^(q)Z  I  X  <q,/2lq,+q/l  >r<rtr(q,|l/2*cl)  f 


q1  li'z  - 


Y  V  ^(k-q’V^kffq;  q'f  I  )  i. 

+  >  > - Li —  ru{q'i  t  ■ 

q'  ij<2  <  ti’+q<ilq’l2>  1  2’ 


q)i 
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(u  -  u£m  (k+q)  +  uj,  (k))  r'  *  (kit;  q  )  = 

=  [nF  (uf  (k+qj  -nr  (uf*  (k»]  <k  +  qf|ki'  >  j-  Uj  (q)  + 

V  <k-q,y|£/(kK,q;  q'lj^  )  _H 

4  <q,  +  qi1|q,i2  > - r  ta' V2 i  «>  ’ 

^  I  2 


(B-5) 


«/  (k)  =  C<<k)-  |<q'l'|kl>|2^(k-q')  (nr(q')  +  «r3  (q'))  (B-6) 

q»f 

To  see  the  connection  with  the  vertex  part  equations,  namely  (3.  3.  6), 
one  first  notices  that  the  right  side  of  (3.  3.  6)  is  independent  of  the  frequency 
part  of  k  ,  the  first  index  of  the  vertex  function, so  that  one  may  assert  that 
the  solution  of  this  equation  is  of  the  form 

r.'(kff';q)  =r.'(k  If';  q)  . 

Then  the  «  integrations  can  be  carried  out  in  the  right-hand  side  and  assuming 
the  G's  to  be  diagonal,  one  gets,  writing  it  in  the  component  form: 
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*  V2 


Now  let 


nF(^(q+q'))  -nF(u^c(q')) 

u-  u/  (q+q1)  +  u/  (q») 

1  f2 


,  r  ,  fu  „  U-  W<  <k+q)  +  «*  (k)  \  r.°°  (k II';  q) 

1  _  (nF(u;(ktq))-„F  ic/'o,,,  )  JWTi-na'?-  'B-8> 

(u  -«/  (k+q)  +  u/  (k))  r.0'0'  (  kil'  ;  q)  = 

=  [nF(u/<r  (k+q))-nF(u)^  (k))  ]  <  k  +  q  i |  kf  >  ^  (r.)^  - 

-  6aa^(q)  £  £  ^(<q'l2|q'fqf1  >)  I?"{q'  ^  ^  ;  q  )  + 


(Ii,(a,ilf2»  <*»<„,*  (B-7) 
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2/(k-q')-2^{kil'q;  q'/j^) 
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which  are  identical  in  structure  with  those  derived  above  (B  -  3, 4,  5)  except 
for  the  external  field  terms  which  have  been  swallowed  in  the  Green's  function 
technique  as  the  vertex  part  involves  a  variational  derivative  with  respect  to 
these  aid,  hence,  their  presence  here  appears  as  (The  notation  used 

here  is  very  slightly  different  from  the  text  and  should  cause  no  great 
confusion .  ) 


T '  r- 
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APPENDIX  C 


OVERHAUSER  SOLUTIONS  INVOLVING  BLOCH  ELECTRONS 

-  f  ■  ■  m3 EH  ^ 

In  the  Epilogue,  the  fifth  problem  mentioned  concerns  the  extension 

of  the  results  of  Overhauser  [18j  obtained  in  Section  II- 3  to  Bloch  electrons. 

This  is  particularly  of  use  since  the  Chromium-Rhenium  system  seems  to  show 
anSDW  character  demanding  a  two-band  model.  This  was  brought  to  the 
attention  of  one  of  us  {A.  K.  Rajagopal)  by  Dr.  Van  Zandt  of  Lincoln  Laboratories. 
In  this  appendix  we  give  a  generalization  of  our  results  of  the  Stoner  theory  for 
Bloch  electrons  given  in  Section  II — 4  so  as  to  give  an  Overhauser-typc  theory 
for  Bloch  electrons.  The  results  of  Section  II- 4 follow  when  Q  =  0  and  those  of 
Section  If  -  3  follow  when  the  Bloch  functions  are  replaced  by  plane  waves,  and 
the  band  indices  omitted.  The  results  obtained  here  are  all  formal. 

In  order  to  arrive  at  a  suitable  generalization,  we  assume  that  the 
electron  of  spin  up  and  wave  vector  k  in  band  l  is  associated  only  with  the 
one  of  spin  down  in  the  same  band  but  with  a  vector  (k  +  Q)  .  With  this 
prescription,  the  following  redefinitions  of  the  expression  (2.  3.  1)  are  made: 


■I  S 


V  IBZ 


G.  1(11') 


G, |(11') 


■is 

f  IBZ 

•i  i 


Gi.(ll') 


V  IBZ 


■If 


I'  IBZ 


bf'k'  Gj  j  (i,k’;ti  -  tp) 


hf'k'+Q*1*  hf'k'+Q*1'*  G|  1  ^'k’>  1 1 -t i * > 


h/'k’*1*  hf'k'+Q*1*  Gf|  (^’k*  i  4i  -tp) 


’bf'k'+Q^  hf'k'^1')  G|  (U'k'j  -tii) 


(C-l) 
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Then  a  calculation  similar  to  the  one  given  in  Section  II- 3  follows  and  we  giv# 


here  only  the  final  results. 

vT 


clt(U)  =  ei(K)  +  l£  fatZ  C„VlCiU').l*+ 


V  IBZ 


c(,  (k)  =  <,<k+QMi£ 

V  IBZ 


V  IBZ 


«I,W‘X  K+°£ia  Glt<(i'k'>u',e 


iu'o+  d3k'du' 


V  IBZ 


Here  e^(k)  is  the  energy  of  the  electron  in  band  f  ,  and  the  IBZ  indicates 
that  the  integration  in  k  space  is  confined  to  the  first  Brillouin  zone.  As 
in  Section  II - 4, 


■  i  S’  \k,(‘>  v2(2)V<l-2>  v3|2)  V4,i>d3id3z 


3.  .3. 


Moreove  r, 


Ul-w  =  r«°r^) +  +eu(/k)git(fk> 


(C-3) 


i 
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Thu«#  in  Section  II- 3  we  may  replace*  Oy  by  so  that  one  has  finally 

G  (Ik » 


fcoa2^ w)] 


coa^ainfl^jjj^fg^fk;  w^gjfkju) 


[co.^.^/m^-  fg+^i  u)-g_(fk;  cj)]  [•in20&g+(fk;  uj+coa^^gflk;  u)] 


whore 


g+  {/k;  u)  =  1/tu  -  U|-  (k)]  . 


(C-4) 


(C-5) 


Finally  then 


ei|(k)  =  ey(k)-^  y  [co82ai(k,nF(u/{k»))+.in20i,klnr(uir(k'))) 


f  IBZ 


U\M  -  ',<k+Q>  -I  J  ■^^SSrkH2t>i"29rk’"F(“/t'<k'»tc':'“2e/.k' 

^  ibz  1  ' 

V  P  d^k'  tkl'k'  +  Q  /g ♦  I  )  , 

gtJ^  2/ J  7TI73- ^fk+Q /'k'  N/gTTrrFy  COB°i'k'81n0rkl{nF(tJi'(k,))'nF((Ji1  {k,))J 


i'  IBZ 


?if(£k>  =  ^  I*  .3  ^  Ik/'k'+Q  / gjj(f'k')  CO80/«K'8in5rk^rF^WJ^k,^'nF^u,l^k'^^ 


i'  IBZ 
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l.l  -  A 


The  equations  for  gjj,  gj  |  may  be  written  more  expHcitlynoijosfi  ni  ^4riV 


-r  t*\  T  f  d3k'  V^'fc'+Q  g|l(i'k,)  rn  (M+Ikru^lu'lk* 

^  8tl(£k)-lJ  “T^fk+Qi-k-  ,  +/k>)  -u-  kT]  F  f  F 

/  *  •*  .»*  •  'i*  i 

-  * 1  “  f  i  *  r-  *  , 

C  and  a  similar  equation  for  g|  j  (£k).  , 

[  X 

Lastly, the  expression  for  total  energy  of  this  system  in  HF  may  be 
j  f  'v  )fi.‘  1  • 

computed  using  these  Green's  functions  as  before.  This  is  seen  to  be  » ' 


(**))] 

7) 


;}  *)  V  . ' 


—  ^  j*  -^3  [(e/qjcos^*  e^q+Osin^nytu/tq)) 
fio  t  IBZ  W 

+  (e^q^in^^+e^q+QJcoa^^)  nF(Uj~(q))l 

-ini'  ?S>“Vr<-Iw>* 

IV  T^2.  14**  L  /  ,  , . 

+  sin  (q)))(cos20FklnF(u^(q'))+3in20Fk,nF(u>F  ( q* >)> 

^i'q'  +  Qfq+Q  (3in\nF(^))tcosVrU<(q)))(8in2<,<'k'nF,Ui’  (q,))  + 

+  co320(,k,  nF(ut,  (q'))) 

+‘2//,q'i^Q  j!"!-1  cosO  sin 0 .  cos  9  tsin0  , 

^i'q'rQ  Iq  sj  gj  |  { fq)  g|  |  «' q' )  *q  ^  *'q'  *V 

[nF(wi+(q))-nF(u‘  (q))][nF(oiJ'(q'))-nF(u)i,'(q'))] 


+  /?yrq'+Q/q  /8ltUq)gn(i,q,) 

cos  0.  sin  0.  cos  8  r  ,  sin  0.,  , 
,q  q  q  q 

n 

[nF(Ui+('4))-rlF{u/'!q!)]  [nf (<J^(q'))-nF(<jF(q'))] 
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This  generalizes  <2.  5.  26)  for  Bloch  electrons.  We  may  remark  that  just  as 
for  Q  =  0,  the  equations  for  g|  j  and  gj  j  are  now  coupled  integral  equations. 

In  the  above  calculation,  the  Umklapp  processes  were  neglected  which  implies 
the  neglect  of  local  field  corrections  and  thus  broad  bands  [44).  This  formally 
completes  the  analysis  of  SDW  for  Bloch  electrons. 
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APPENDIX  D 


P-STATE  INSTABILITY 


The  solution  obtained  in  Section  III-2  was  not  quite  suitable  for  discussing 
the  P  state  stability  when  the  Coulomb  case  for  f  *  0  is  considered.  Here  we 
derive  the  correct  criterion  for  P  state  instability.  Consider  the  spin  wave 
equation  (3.  2.42)  in  the  limit  u  =  0  f  5  =  0  (P  state),  lhen 


(l  ~r2'v,p-k)[ni 


.{u(k+q))-n.,  (u  (k))] 


>  T(p;q) 


■j 


2/y(p-k)[nF{u(k+q))-nF(w{k))]  T(k;q) 


(D-l) 


Now  let 


r (p;q)  =  ^  ql  Y  tm 


i,  f,m 


2X.  2  ,  _2  .  2 

^  p  +  k  +  £  k  _  a  «iA 

^y(p'k)  =  W  °f(  - ZpE - }  Yfm(p)  Yfm(k) 

tm 


(D-2) 


nf  (u  (k  +  q))- nF  (w  (k))  =  ^  q  '  'i 

;  tm 


'«  "ta  (k» 


in  (D-l)  and  integrate  over  p  .  Then  we  get 
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i2iim  \ 


Equating  the  coefficients  of  the  like  powers  of  q  we  get  since 

n  <°>  =  0 
fm 

2  L  .  /  p2+k2  +  e2k  2 

-£<%-)' /l  r,W(p)t  fjr  X  IkdkQ<  — 5s 
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fm 


It  is  easy  to  prove  that 

k  2 

n,  (1)(k)  =.  -  <^)1/2  lv(k)'6(-^r-u(k))fi^ 


‘in  'T'  '  k 

Then  (D  5)  takes  the  form 


mo 
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Let  us  take  u  ~  k*  /2m  so  that  on€  finally  gets  a  condition  for  nontrivial  solution 
of  (D-7)  as 

ar 


i  =  -p5  [Qoo  +  r/2)  -  Qj  {i  +  r/2>  ] 


(D-8) 


,Now, 


Qq(x)  -  Qj  (x)  = (1-x)  In  |  +1 


or 


or 


ar  . 

i  =  ri  -  ^  (i  +  )  ] 

Q  j*  2 

1  =  [1  -  Jp-fn  (1  +  )  ] 


(D-9) 


Note  that  this  is  precisely  (2.  5.  16)  when  evaluated  for  f  =  0.  For  £  -*  0  ,  this 

gives  the  Bloch  condition  ar  =  V  .  For  f  -*  *",  fn  (1  +  ~  — ly-  -  — ? 

—  s  ^  | 


and 


ar 


[i  -  i  +  -Vi  --  ( 


ar 


)  2  , 


Z  3  *  2 

but  a rs/jt  ^  =  y  ( - )  ,  so  that  wc  get  the  Stone  r  c  rite  r ion  K O' /ep  =  y 

€F 
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